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ABSTRACT 


Thermophoto  voltaic  (TPV)  energy  conversion  utilizes  photons  from  a  thermal 
radiator  to  convert  photonic  energy  to  electrical  energy.  Due  to  the  nature  of  the  system,  the 
thermal  radiator  must  emit  uniform  radiation  and  therefore  maintain  a  uniform  temperature 
profile  in  order  to  achieve  maximum  efficiency.  Heat  pipe  technology  can  effectively  meet 
the  demand  for  an  isothermal  emitter  as  it  utilizes  near  isobaric  phase  changes  to  transfer  heat 
at  a  uniform  temperature.  In  this  project,  heat  pipes  are  explored  for  use  in  TPV  energy 
conversion  systems.  A  flat  heat  pipe  offers  many  advantages  over  the  conventional 
cylindrical  design.  These  include  increased  surface  area  to  volume  ratio  in  order  to 
maximize  power  density,  as  well  as  the  ability  to  stack  or  layer  the  system  with  photovoltaic 
(PV)  cells  on  both  sides  of  the  flat  heat  pipe  to  utilize  available  energy.  Not  only  do  flat  heat 
pipes  present  unique  engineering  demands  inherent  in  their  operation  to  counteract  pressure 
differences  across  their  vessel  walls,  but  they  are  also  difficult  to  construct.  To  date,  only 
limited  analyses  of  their  thermal  characteristics  have  been  done  for  use  in  performance 
predictions.  Therefore,  it  is  necessary  to  conduct  analyses  to  enable  consideration  of  heat 
pipes  for  implementation  into  TPV  systems.  This  report  details  the  design  and  construction 
of  a  flat  heat  pipe  analysed  both  in  symmetric  and  asymmetric  heating  conditions,  involving 
a  low  temperature  version  of  future  emitter  designs  due  to  safety  considerations.  Water  was 
used  as  a  working  fluid  instead  of  the  liquid  metal  required  to  achieve  the  temperatures  of 
a  functional  emitter.  Despite  this  difference  in  working  fluid,  the  data  presented  is  valuable 
to  both  TPV  and  heat  pipe  research. 

KEYWORDS:  heat  pipe,  emitter,  isothermal,  thermophotovoltaic 
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General  Terms 

area  (usually  cross-sectional)  (m^) 

acceleration  (m/s^) 

constant  value  in  an  equation 

integration  constant 

function  of,  or  friction  coefficient 

function  of,  or  friction 

acceleration  due  to  gravity  (m/s^) 

Grashof  number 
height  (m) 

heat  transfer  convection  coefficient  [W/(m  K)] 
permeability  (m^) 
thermal  conductivity  [W/(m  K)] 
length  (m) 

length  of  subsection  (m) 
mass  (kg),  term  in  fin  equation 
mass  flow  rate  (kg/s) 

Mach  number 
Screen  mesh  number 
number  of  repeated  elements 
Nusselt  number 
pressure  (N/m^),  perimeter  (m) 

Prandtl  number 
heat  rate  (W) 

gas  constant  for  given  vapor  [kJ/(kg  K)] 
radius  (m),  radius  of  curvature  (m) 

Rayleigh  number 

Reynolds  number 

spacer  distance  (m) 

temperature  (K) 

time  (s),  thickness  (m) 

overall  heat  transfer  coefficient  [W/(m^  K)] 

volume  (m^),  velocity  (m/s) 

specific  volume  (m^/kg) 

width  (m) 

Weber  number 
cartesian  coordinate 
cartesian  coordinate 
cartesian  coordinate 
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A, 

As 

A.. 

A. 

C, 

C2 

Co 

Dv 

E, 

Ebx 

F, 

Fv 

g 

hf. 

K' 

k| 

Fefr 

ks 

kv 

k^ 

Le 

FefT 
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m, 
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Po 

Qb 

Q.n. 

Qs 

Q.S 

QL 

<ie 

P-h.w 

Rv 
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T 

amb 

Th 

Flocal 


Specific  Terms 

cross-sectional  area  of  liquid  flow  path  (m‘) 
total  wetted  surface  area  of  wick  (m^) 
vapor  flow  cross-sectional  area  (m^) 
area  in  x-direction  (m^) 

first  radiation  constant  (3.7415  x  10''®  W  m^),  first  constant  of  integration 
second  radiation  constant  (1.4388  x  10"^  m  K),  second  constant  of  integration 
speed  of  sound  evaluated  at  evaporator  end  cap  temperature  (m/s) 
diameter  of  vapor  core  (m) 
total  emissive  power  (W/m^) 

blackbody  emissive  power  at  absolute  temperature  T  for  a  single  wavelength 
(W/m^) 

liquid  frictional  coefficient 
frictional  coefficient  for  vapor 
acceleration  due  to  gravity  (m/s“) 
latent  heat  of  vaporization  (kJ/kg) 
ratio  of  specific  heats 

thermal  conductivity  of  working  fluid  in  liquid  phase  [W/(m  K)] 
effective  thermal  conductivity  of  the  wick  [W/(m  K)] 
thermal  conductivity  of  wick  material  [W/(m  K)] 
thermal  conductivity  of  vapor  [W/(m  K)] 
thermal  conductivity  of  vessel  wall  material  [W/(m  K)] 

Evaporator  length  (m) 

effective  fluid  return  length  (m) 

total  length  (m) 

liquid  mass  flow  rate  (kg/s) 

vapor  mass  flow  rate  (kg/s) 

pressure  at  evaporator  end  (N/m^) 

maximum  heat  rate  at  boiling  limit  (W) 

maximum  heat  rate  at  entrainment  limit  (W) 

maximum  heat  rate  at  sonic  limit  (W) 

maximum  heat  rate  at  viscous  limit  (W) 

heat  transport  capability  (W  m) 

maximum  heat  flux  into  evaporator  section  (W/m^) 

hydraulic  radius  of  wick  surface  pore  (m) 
inner  radius  of  pipe  wall  (m) 
vapor  space  radius  (m) 

liquid  velocity  through  wick  screen  pores  (m/s) 
temperature  of  ambient  environment  (K) 
temperature  of  hot  (radiating)  body  (K) 
local  temperature  (K) 
steady-state  temperature  (K) 
vessel  wall  thickness  (m) 
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wick  thickness  (m) 

^Pcap.max  maximum  capillary  pumping  pressure  (N/m’) 
Ap^  pressure  drop  due  to  gravitational  force  (N/m') 

Api  liquid  pressure  drop  (N/m^) 

Ap^.  pressure  drop  in  vapor  (N/m^) 

AT^nt  Wall  Superheat  temperature  difference  (K) 


Greek  Symbols 

6^.  thickness  of  vapor  flow  cross-section  (assume  uniform  for  entire  heat  pipe)  (m) 

P  coefficient  of  thermal  expansion  (1/K) 

e  surface  emissivity 

X  wavelength  (m) 

wavelength  at  which  an  ideal  radiator  emits  maximum  amount  of  energy  , 

Pi  liquid  viscosity(N  s/m^) 

Po  vapor  density  evaluated  at  evaporator  end  cap  temperature  (kg/m^) 

Pi  density  of  liquid  (kg/m^) 

p^,  vapor  density  (kg/m^) 

o  liquid  surface  tension  (N/m),  Stefan-Boltzmann  constant  [5.67  x  10'*  W/(m^  K"*)] 

4)  angle  of  heat  pipe  from  horizontal  (degrees) 

(p  wick  porosity 
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1.0  INTRODUCTION 

This  project  report  entails  the  design,  construction,  and  both  the  theoretical  and 
experimental  thermal  performance  investigation  of  a  flat  heat  pipe  using  water  as  the 
working  fluid.  The  project’s  findings  will  help  determine  the  viability  of  implementing  flat 
heat  pipes  into  thermophotovoltaic  energy  conversion  systems. 

1.1  OBJECTIVES 

The  flat  heat  pipe  was  designed  to  meet  the  following  criteria  (as  determined  by  both 
the  sponsors  of  the  project  and  the  research  team); 

1.  Design  and  construct  the  heat  pipe  entirely  out  of  Monel  400. 

2.  Exhibit  two-dimensional  heat  transfer  behavior. 

3.  Feature  ultimate  dimensions  of  4  feet  long  (1.22  m)  by  1  foot  wide  (30.48  cm)  by 
inch  thick  (1.27  cm). 

4.  Maintain  an  operating  temperature  of  up  to  212°F  (100°C). 

5 .  Possess  a  means  of  internal  support  for  structural  rigidity  without  sacrificing  two- 
dimensional  flow  capabilities. 

6.  Operate  without  the  presence  of  any  noncondensable  gases,  such  as  air,  inside  of 
the  heat  pipe  (gas  loading). 

The  following  experimental  and  theoretical  goals  were  set: 

1 .  Monitor  the  condenser  region  temperature  profile  as  closely  as  possible,  via 
thermocouples  and  infrared  video  camera. 

2.  Test  asymmetric  heating  conditions  and  different  heat  inputs. 

3.  Test  performance  at  different  heat  inputs,  fluid  charges,  and  varying  inclination 
angles. 

4.  Modify  important  heat  pipe  equations  for  flat  geometry. 

5.  Model  experimental  and  theoretical  heat  pipe  performance  against  a  conventional 
fin  of  the  same  dimensions. 

6.  Analyze  experimental  data  to  determine  overall  performance  and  feasibility  of 
flat  heat  pipes  for  a  TPV  energy  conversion  system. 

1.2  METHODOLOGY 

The  design,  experimental,  and  theoretical  goals  were  met  in  the  following  manner: 
Step  1 :  Research  and  Design 

A  literature  search  was  conducted  to  determine  the  mechanisms  of  heat  transfer  in 
heat  pipes  and  how  they  might  differ  in  a  flat  heat  pipe  from  a  cylindrical  heat  pipe.  Several 
experts  in  the  field,  with  specialties  ranging  from  design  and  theoretical  analysis  to 
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fabrication,  were  consulted  and  their  recommendations  considered.  Suppliers  of  the 
necessary  materials  were  located.  The  final  design  was  completed  incorporating  all  of  the 
above  research. 

Step  2:  Construction  and  Testing 

A  subscale  model  of  the  final  design  was  submitted  to  Technical  Support  Department 
(TSD)  at  USNA  for  construction  to  allow  fabrication  difficulties  to  be  solved  before 
attempting  the  final  design.  Instrument  calibration  and  basic  heat  pipe  operational 
procedures,  such  as  fluid  charging,  were  perfected  with  this  model.  The  subscale  model  also 
served  to  provide  insight  that  was  incorporated  in  revisions  for  the  fullscale  model.  Finally, 
the  fullscale  heat  pipe  was  constructed,  tested,  and  analyzed  under  the  test  conditions. 

Step  3:  Analysis 

The  raw  data  from  the  experiments  were  collected  and  reduced.  The  infrared  camera 
data  was  used  to  aid  in  explaining  the  purely  quantitative  thermocouple  data  as  well  as 
provide  visual  understanding  of  heat  pipe  phenomena.  The  temperatures  of  regions  at  the 
same  distance  axially  from  the  heater  end  were  compared  to  other  equidistant  regions  and 
analyzed  for  thermal  performance  assessment  at  varying  heat  inputs,  fluid  charge,  and  angles 
of  inclination.  A  traditional  fin,  with  the  same  dimensions  as  the  fullscale  heat  pipe,  was 
modelled  against  the  heat  pipe’s  performance.  This  demonstrated  the  heat  pipe’s  ability  to 
transfer  large  quantities  of  heat  with  a  near  isothermal  profile;  a  quality  not  found  in  a 
traditional  fin. 
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2.0  THERMOPHOTOVOLTAICS  (TPV)  BACKGROUND 

TPV  energy  conversion  is  a  phenomenon  by  which  infrared  thermal  radiation  is 
converted  into  electrical  power  with  no  moving  parts.  The  basic  components  of  the  system 
include  the  thermal  radiator  and  a  semiconductor  diode  (thermophotovoltaic  cell)  that 
converts  thermal  radiation  into  electrical  power.  Only  photons  with  energy  greater  than  the 
bandgap,  or  threshold  energy  of  the  TPV  cell,  can  be  converted  into  electrical  energy.  The 
radiation  not  converted  either  creates  heat  in  the  cell  or  is  recycled  to  the  radiator  via 
reflection  if  such  a  system  (known  as  spectral  control)  exists.  The  bandgap  energy  is 
determined  by  the  TPV  cell  material  and  is  not  a  variable  parameter  (Erickson,  1997).  For 
a  schematic  of  the  basic  TPV  system,  see  Figure  1. 


Spectral  Control 


TPV  Cell 


Figure  1.  Basic  Thermophotovoltaic  Energy  Conversion  System 


Figure  2  shows  the  emissive  power  over  differing  wavelengths  for  theoretical 
blackbodies  at  different  temperatures.  A  blackbody ,  or  ideal  radiator,  is  an  object  that  emits 
the  maximum  possible  amount  of  radiation  at  any  given  wavelength.  The  emissive  power 
is  a  function  of  only  the  temperature  of  the  emitter.  The  emissive  power  for  a  single 
wavelength  is  given  by  Planck’s  law  in  equation  2.1: 


Cl 


(2.1) 


From  Figure  2  it  is  apparent  that  increasing  temperature  will  increase  the  emissive 
power  of  the  blackbody  for  a  given  wavelength.  The  wavelength  at  which  the  emissive 
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power  is  maximized  (which  corresponds  to  the  peak  of  the  curve  for  each  temperature  in 
Figure  2)  is  found  by  Wien’s  displacement  law  given  in  equation  2.2; 


(2.2) 


,  _  2.898x1 K 

'l-max“  7’ 


Monochromatic  Emissive  Power 


Figure  2.  Planck’s  Law 


The  total  emissive  power  of  a  blackbody  at  a  given  temperature  is  the  area  underneath 
the  curve  corresponding  to  the  given  temperature  in  Figure  2.  This  is  represented  in  equation 
2.3,  the  Stefan-Boltzmann  law  (where  o  is  the  Stefan-Boltzmann  constant): 


E,(T)  =  oT^ 


(2.3) 


Since  power  is  proportional  to  temperature  to  the  fourth  power,  small  changes  in  the 
temperature  of  an  emitter  can  have  great  effects  on  emissive  power. 


2.1  PHOTOVOLTAIC  PHYSICS 

The  physical  mechanism  within  the  TPV  cells  that  governs  the  conversion  of  photons 
into  electrical  energy  is  the  photoelectric  effect,  the  same  principle  which  standard 
photovoltaic  (PV)  cells  (also  known  as  solar  cells)  utilize.  The  active  layer  of  these  cells  is 
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constructed  of  a  semiconductor  material  with  a  P-N  (Positive-Negative)  junction.  The 
material  used  to  produce  the  active  cell  layers  determines  the  bandgap  energy,  or  the  energy 
necessary  to  dislodge  an  electron  from  the  material’s  outermost  electron  shell.  As  photons 
from  the  thermal  radiator  travel  into  the  TPV  cell,  they  travel  through  the  “window”  and  are 
absorbed  by  the  semiconductor.  If  the  energy  of  the  photon  striking  the  semiconductor  is 
sufficient,  the  semiconductor  will  eject  an  electron  from  its  outer  shell,  leaving  a  “hole.” 
If  not,  it  will  cause  the  temperature  of  the  material  to  rise.  If  the  difference  in  energy 
between  the  “hole”  and  the  ejected  electron  is  greater  than  the  bandgap  energy  of  the 
material,  the  pair  (electron  and  hole)  will  cross  the  voltage  barrier  of  the  cell.  When  this 
charge-carrying  pair  moves  towards  similarly  charged  regions  and  moves  away  from  layers 
of  unlike  charge,  it  becomes  an  electrical  current,  or  usable  electrical  energy.  The  amount 
of  energy  available  to  be  converted  into  electricity  by  the  photoelectric  effect  is  determined 
by  Planck’s  law  (see  equation  2.1).  Only  photons  that  have  energy  (emissive  power)  greater 
than  bandgap  energy  can  produce  an  electron-hole  pair  (and  therefore  produce  electricity) 
(Courts  and  Fitzgerald,  1998). 

As  established,  PV  (and  therefore  TPV)  cells  can  only  convert  certain  wavelengths 
of  electromagnetic  radiation  into  usable  energy.  Therefore,  the  system  efficiency  is  a 
function  of  the  wavelengths  emitted  by  the  thermal  radiator.  As  shown  in  Figure  2, 
wavelength  is  a  function  of  temperature.  Therefore,  TPV  cell  efficiency  is  a  function  of 
radiator  temperature.  Methods  do  exist  of  reflecting  back  emitted  wavelengths  that  cannot 
be  utilized  by  the  TPV  cells  for  reabsorption  into  the  emitter.  This  is  known  as  spectral 
control.  It  ensures  that  the  TPV  cells  receive  only  usable  wavelengths  of  photons,  and  results 
in  higher  energy  conversion  efficiencies.  However,  spectral  controls  cannot  moderate  the 
amounts  of  radiation  that  come  from  different  regions  of  the  emitter  (corresponding  to 
locations  on  the  emitter  of  different  temperature)  (Borowsky  and  Dziendziel,  1996). 

2.2  TPV  CELL  EFFICIENCY 

TPV  and  PV  cells  must  be  linked  together  in  series  to  produce  an  appreciable  voltage, 
as  the  voltage  output  of  a  single  cell  is  small.  Since  the  cells  are  in  series,  the  same  amount 
of  current  must  flow  through  each  cell.  Experience  has  shown  that  if  one  cell  in  a  series  of 
TPV  cells  receives  fewer  photons  above  bandgap  energy  than  the  other  cells  (due  to  a  low 
temperature  region  on  the  emitter)  then  the  performance  of  the  entire  series  of  cells  is 
downgraded  (Lindler,  1999).  Therefore,  the  emitter  must  exist  at  a  uniform  temperature  for 
maximum  efficiency.  This  presents  a  substantial  problem  for  emitter  design  in  TPV  energy 
conversion  systems.  Combustion  processes  can  have  inherent  non-isothermal  profiles  and 
therefore  can  not  operate  TPV  cells  at  a  high  efficiency.  However,  heat  pipes  provide  a 
means  of  heat  transfer  that  maintains  a  nearly  isothermal  surface  and  could  serve  as  emitters 
in  TPV  generators.  Figure  3  shows  just  such  a  system.  Heat  pipes  are  heated  axially  via 
fossil  fuel  combustion.  Alternating  heat  pipes  and  TPV  cells  results  in  a  large  surface  area 
to  volume  ratio. 
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Combustion  Chamber 

Figure  3.  USNA  TPV  axial  flow  concept 
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3.0  HEAT  PIPE  INTRODUCTION 

Heat  pipes  have  existed  for  decades  and  stem  from  a  simple  principle:  vapor 
evaporates  and  condenses  at  a  uniform  temperature  for  a  given  pressure.  A  heat  pipe  consists 
of  a  hermetically  sealed  vessel,  a  working  fluid,  and  a  wick  that  lines  the  vessel.  The 
working  fluid  is  vaporized  (by  an  external  heat  source)  and  leaves  the  wick  as  a  vapor.  The 
vapor  travels  from  the  evaporator  end  to  the  condenser  end  of  the  heat  pipe  and  condenses 
in  the  wick,  releasing  heat.  Due  to  the  change  in  phase  of  the  working  fluid,  the  heat  is 
delivered  to  the  condenser  end  at  an  almost  isothermal  profile.  The  fluid  returns  to  the 
evaporator  end  from  the  condenser  end  via  the  capillary  pumping  action  of  the  wick.  Figure 
4  shows  the  basic  configuration  of  a  cylindrical  heat  pipe. 


Evaporator  End  Adiabatic  Region  Condenser  End 


Figure  4.  Basic  Heat  Pipe  Configuration 


Constructing  and  operating  a  heat  pipe  is  an  involved  process,  with  many  limitations 
to  be  considered.  First,  the  fluid  and  vessel  must  be  of  high  purity.  If  the  fluid  is  not  pure, 
then  its  heat  transfer  properties  may  fluctuate  from  the  expected  standard.  Impure  fluid  can 
also  deposit  undesirable  substances  in  the  heat  pipe  vessel  and  vapor  space,  such  as 
noncondensable  gases.  Undesired  gases  in  the  vapor  space  serve  to  effectively  blanket  the 
condenser  area  and  raise  the  temperature  of  the  heat  pipe  in  the  areas  removed  from  the  gas 
while  the  areas  where  the  gas  is  located  remain  colder  than  the  rest  of  the  heat  pipe.  This 
effect,  called  gas  loading,  differs  for  varying  heat  inputs.  Figure  5  serves  to  illustrate  this 
effect  (Ashcroft,  1996). 

If  the  vessel  is  impure,  then  the  working  fluid  may  possibly  erode  parts  of  the  vessel 
(e.g.  impurities)  (Ashcroft,  1996).  The  vessel  material  must  also  handle  large  compressive 
and  tensile  forces  (Peterson,  1994). 
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'O  Non-functioning  condenser  area 

Low  Power 

Medium  Power 

High  Power 

Evaporator  Condenser 

Figure  5.  Gas  Loading  Effect 

3.1  STEADY-STATE  LIMITS 

The  heat  pipe  has  many  different  operating  limits  for  given  input  heating  powers  and 
operating  temperatures.  When  too  much  heat  is  applied  the  working  fluid  is  superheated  and 
bubble  formation  occurs  in  the  wick,  which  then  dries  out  and  interrupts  heat  transfer.  This 
is  known  as  the  boiling  limit  and  is  represented  by  equation  3.1  (for  cylindrical  heat  pipes) 
(Faghri,  1995): 


(3.1) 


The  viscous  limit  primarily  exists  in  heat  pipes  that  have  startup  conditions  in 
which  the  working  fluid  is  solidified  (such  as  in  high  temperature  liquid  metal  heat  pipes). 
It  is  possible  for  the  working  fluid  to  liquefy  and  be  vaporized  from  the  evaporator  region 
and  condense  in  the  condenser  region  without  sufficiently  heating  the  frozen  fluid  in  the 
condenser  region  to  allow  flow  back  into  the  evaporator  region.  The  evaporator  will  then  dry 
out  and  the  heat  pipe  will  not  reach  normal  operating  conditions.  In  general,  the  viscous 
limit  can  be  reached  when  a  heat  pipe  is  operating  at  temperatures  below  its  normal  operating 
range,  when  viscous  forces  dominate  the  vapor  flow  and  vapor  pressure  in  the  condenser  end 
is  reduced  to  zero,  limiting  the  heat  transport  capabilities  of  the  heat  pipe.  The  viscous  limit 
is  expressed  by  equation  3.2  (  for  cylindrical  heat  pipes)  (Peterson,  1994): 


64//^ 

At  low  vapor  densities  (and  hence,  low  temperatures),  the  vapor  velocity  must  be 
high  to  transfer  the  necessary  heat.  For  a  constant  area  duct,  the  vapor  velocity  can  not 
exceed  the  local  sonic  velocity.  Therefore,  problems  such  as  dryout,  shock,  and  large 
temperature  fluctuations  result.  The  sonic  limit  is  related  as  follows  in  equation  3.3  (Dunn 
andReay,  1976): 


Po^o^/g^v 

yl2{K'+l) 


(3.3) 


In  certain  wick  designs,  it  is  also  possible  for  the  vapor  that  is  traveling  to  the 
condenser  end  to  capture,  or  entrain,  the  returning  liquid  and  cause  dryout.  This  limit, 
known  as  the  entrainment  limit,  is  expressed  as  (Chi,  1976): 


(3.4) 

The  most  important  limit  for  design  considerations,  the  capillary  limit,  is  met  when 
the  vapor  and  liquid  pressure  drops  exceed  the  maximum  capillary  head  that  the  wick 
produces.  This  maximum  head  is  related  to  the  wick’s  effective  pore  size.  The  capillary 
limit  occurs  when  the  pumping  ability  of  a  given  wicking  material  to  provide  the  circulation 
for  the  working  medium  can  no  longer  return  fluid  as  fast  as  the  evaporator  vaporizes  it. 
This  limit  will  cause  the  evaporator  end  to  dry  out,  resulting  in  damage  to  the  vessel  and 
wick.  The  capillary  limit  is  the  highest  heat  rate  that  most  heat  pipes  can  operate  at  to  achieve 
minimum  design  operating  temperature  and  is  therefore  usually  the  only  limit  given  for  heat 
pipes.  See  Figure  6  for  the  relations  of  temperature  and  heat  rate  to  heating  limits  for  a 
conventional  cylindrical  heat  pipe.  Equation  3.5,  the  equation  for  the  capillary  limit  for  a 
cylindrical  heat  pipe,  neglecting  the  pressure  drops  due  to  evaporation  and  condensation,  is 
(Faghri,  1995): 


(3.5) 
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(QL)^3p,„ax  is  the  heat  transport  capability  (W  m),  F,  is  the  liquid  friction  loss  coefficient,  and 


Fy  is  the  vapor  friction  loss  coefficient.  Equations  3.6  through  3.10  further  define  the 
elements  of  the  capillary  limit  equation  (equation  3.5). 


^Pcap,max~  (2o)/rjff 


Apcap,max  the  maximum  capillary  pumping  head 


APcap,max  ^  Ap,  +  Ap^  +  Apg  where 


^Pl  = 


Ap,  is  the  pressure  drop  in  liquid 


Pj'effrh^ 

pXA 


Apy  is  the  pressure  drop  in  the  vapor 


APg  =  p,  gL,  sin  ({) 


(3.10) 


Ap  is  the  pressure  drop  due  to  gravity  and  sin  (})  =  angle  of  inclination  from  horizontal 
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Figure  6.  Heating  Limits 
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Figure  7  depicts  the  wetting  angle  for  a  fluid  and  its  interface.  The  maximum 
capillary  head  is  calculated  under  the  assumption  that  the  working  fluid  is  perfectly  wetting 
the  wick  material.  This  assumption  stipulates  that  the  contact  angle,  0,  between  the  fluid  and 
wick  is  equal  to  zero.  A  contact  angle  of  zero  means  that  the  fluid  meniscus  is  a  perfect 
hemisphere,  which  is  almost  never  the  actual  case,  as  most  environments  possess  a  wetting 
angle  greater  than  zero.  Therefore,  the  actual  maximum  capillary  head  is  less  than  the  value 
obtained  by  the  equation  3.6.  An  accurate  value  for  maximum  capillary  pumping  head  is 
best  obtained  by  experimental  methods  (Peterson,  1994). 


Sketch  of  side 
view  of  meniscus 


Figure  7.  Wetting  angle 


3.2  WICK  DESIGNS 

For  heat  pipes,  there  are  many  different  wicking  structures  that  could  be  used  to 
provide  the  path  for  the  return  of  the  condensate  from  the  condenser  back  to  the  evaporator. 
To  obtain  the  desired  heat  transport  characteristics,  the  wick  design  must  be  optimized.  This 
involves  a  trade-off  of  various  wick  characteristics.  To  achieve  a  large  capillary  head,  the 
wick  must  have  a  small  pore  size,  as  shown  in  equation  3.6.  However,  in  order  to  minimize 
the  liquid  pressure  drop  from  the  condenser  end  to  the  evaporator  end  the  wick  should  have 
a  large  pore  size.  This  minimizes  F„  in  equation  3.5.  The  effective  thermal  conductivity  of 
wicking  structures  is  also  an  important  consideration.  Due  to  all  of  these  factors  of  wick 
design,  wick  optimization  is  an  involved  process  (Cerza,  1998). 


Since  the  vapor  condenses  on  the  wick  structure,  the  return  liquid  velocity  is 
relatively  low  and  heat  is  transported  through  the  wick  and  vessel  wall  to  the  environment. 
The  effective  thermal  conductivity  is  dependent  on  both  the  working  fluid,  wick  material. 
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and  number  of  pores  per  unit  area.  It  is  the  effective  thermal  conductivity  that  allows  for  the 
calculation  of  the  maximum  wick  thickness  for  desired  operating  temperatures  since  the 
effective  thermal  conductivity  determines  the  temperature  drop  across  the  wick  for  given 
conditions.  Effective  thermal  conductivity  is  found  using  equation  3.11  (Faghri,  1995): 

(3.11) 

^  ki[{ki  +  k;)-{\-(p){kj-k;)\ 


where  porosity,  (p,  is  wick  open  area  divided  by  total  wick  area. 

The  number  of  wicking  structures  in  existence  vary  widely.  However,  they 
can  be  divided  into  two  basic  categories:  homogeneous  and  composite  wicks.  A 
homogeneous  wick  is  constructed  of  a  single  basic  material.  Common  homogeneous  wicks 
are  the  wrapped  screen  wick,  which  consists  of  one  or  more  layers  or  screen  in  contact  with 
the  vessel  wall;  powdered  metal  wicks,  which  are  metal  powders  sintered  to  the  vessel  wall; 
and  the  grooved  wick,  which  consists  of  grooves  machined  into  the  inner  wall  of  the  heat 
pipe.  Homogeneous  wicks  also  can  contain  arteries,  or  tubes  to  provide  liquid  return  paths 
of  very  low  flow  resistance.  Composite  wicks  are  made  up  of  several  wick  materials 
combined  to  form  a  composite  wick  structure.  Examples  of  composite  wicks  are  overlapping 
screens  of  varying  porosities  or  a  groove  wick  with  a  fine  pore  screen  on  top  of  it.  Many 
different  composite  wicks  have  been  developed  for  specialized  applications  (Faghri,  1995). 

The  working  fluid  of  the  heat  pipe  itself  must  also  be  determined.  For  varying 
temperatures  and  heat  pipe  internal  pressures,  different  working  fluids  will  transfer  unequal 
amounts  of  heat  due  to  dissimilarities  in  surface  tension,  density,  enthalpy  of  vaporization, 
and  viscosity.  These  properties  are  grouped  together  in  the  figure  of  merit,  equation  3.12 
(Chi,  1976): 


Equation  3.12  can  be  used  to  determine  which  fluid  will  have  the  greatest  heat  transport 
characteristics  for  the  desired  operating  temperature.  Table  1  shows  various  properties  and 
figures  of  merit  for  various  different  fluids  over  various  temperatures.  Several  computer 
codes  exist  that  can  predict  some  steady-state  operating  limits  of  a  heat  pipe  if  given  the 
material  and  working  fluid  data  for  conventional,  one-dimensional  heat  pipes. 


Table  1.  Figure  of  Merit  Calculations 
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4.0  FLAT  HEAT  PIPE  PROPERTIES 

Flat  heat  pipes  are  very  similar  to  cylindrical  heat  pipes.  The  only  real  difference 
between  the  two  is  geometrical.  While  this  may  seem  a  minor  difference,  it  presents  many 
challenges  from  a  theoretical,  engineering  standpoint. 

Typically,  heat  pipes  are  used  to  transfer  large  amounts  of  heat  across  a  large  distance 
with  only  slight  energy  loss.  The  cylindrical  design  works  well  to  serve  this  purpose. 
However,  when  designing  an  emitter  for  a  TPV  energy  conversion  system,  it  is  advantageous 
to  have  a  large  surface  area  to  volume  ratio  in  order  to  maximize  power  density  of  the 
system.  A  flat  heat  pipe  was  conceived  for  this  purpose.  Flat  heat  pipes  are  a  new 
technology,  with  few  theoretical  models.  Flat  heat  pipes  are  also  relatively  untested  for 
transferring  heat  from  one  end  to  the  other  (as  opposed  to  side  to  side  transfer).  Due  to  the 
different  surface  area  geometry,  flat  heat  pipes  also  have  different  flow  and  structural  design 
considerations  than  those  of  cylindrical  heat  pipes.  This  chapter  quantifies  these  differences 
and  allows  for  consideration  of  these  dissimilarities  in  heat  pipe  design. 

4.1  FLAT  HEAT  PIPE  HEATING  LIMITS 

Flow  properties  in  cylinders  are  much  different  from  those  in  rectangular  geometries 
such  as  flat  plates  and/or  boxes.  The  flow  of  a  thin  film  over  a  flat  sheet  (such  as  in  the 
liquid  return  loop  of  a  flat  heat  pipe)  is  not  the  same  as  the  flow  of  a  cylindrical 
circumferential  film  (which  corresponds  to  a  cylindrical  heat  pipe),  having  different  fiiction 
factors  and  flow  restrictions.  Also,  vapor  flow  through  a  cylindrical  space  differs  greatly 
from  vapor  flow  through  an  object  with  a  rectangular  cross  section,  having  dissimilar  fiiction 
factors  and  flow  restrictions.  The  flow  property  differences  can  present  major  complications 
for  flat  heat  pipe  design,  as  they  alter  steady-state  operating  limits. 

The  limit  most  affected  is  the  capillary  limit.  This  is  significant  because  the  capillary 
limit  is  usually  the  highest  heating  limit  that  a  heat  pipe  is  designed  to  operate  at  in  steady- 
state  conditions.  The  capillary  limit  involves  the  flow  of  the  liquid  in  the  wick’s  return  path 
and  the  pressure  loss  of  the  vapor  and  liquid  through  their  respective  areas  of  the  heat  pipe. 
Also,  since  the  USNA  flat  heat  pipes  do  not  have  adiabatic  sections,  a  more  accurate 
approximation  of  the  mean  liquid  return  length  in  the  wick  was  needed.  The  derivations  of 
the  various  components  of  the  capillary  limit  for  assumed  laminar  flow  in  a  flat,  rectangular 
heat  pipe  are  shown  in  Appendix  1 .  The  capillary  limit  for  a  flat  heat  pipe,  when  only  one 
side  of  the  evaporator  region  is  heated  is  given  in  equation  4.1.  In  equation  4.1  q^  is 
maximum  heat  flux  into  evaporator  section.  Equation  4.2  defines  K,  the  permeability  for  a 
wire  mesh  wick.  To  achieve  a  large  capillary  limit  value  the  effective  pore  size  of  the 
wire  mesh  in  the  wick,  must  be  minimized  to  achieve  a  large  capillary  pumping  head. 
Decreasing  r^fj-  also  decreases  permeability,  which  can  lower  the  capillary  limit  value. 
Proper  flat  heat  pipe  design  involves  optimization  of  both  of  these  factors  to  obtain  a  high 
capillary  pumping  head  while  still  maintaining  an  acceptable  permeability  value.  Equations 
4.1  through  4.5  were  derived  by  Associate  Professor  Cerza  (1999). 
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The  permeability  of  the  wire  mesh  wick  calculated  by  equation  4.2  can  be  approximated  by 
equation  4.3.  Sample  calculations  were  performed  using  both  equations  and  it  was  found 
that  both  equations  yielded  nearly  equal  values.  Therefore,  equation  4.3  is  preferred  due  to 
its  simplicity. 
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There  are  two  heating  configurations  considered  for  the  fullscale  flat  heat  pipe.  The 
first  involves  the  entire  heater  section  mounted  to  one  side  of  the  evaporator  region.  For  this 
configuration,  the  capillary  limit  is  found  by  equation  4.1.  The  other  arrangement  of  the 
evaporator  section  consists  of  heaters  mounted  to  both  sides  of  the  evaporator  region.  For 
this  heater  setup,  the  capillary  limit  is  found  with  equation  4.4. 
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This  capillary  limit  (two  heated  sides)  will  theoretically  be  higher  than  the  limit  with 
only  one  side  heated,  as  the  two  sides  of  the  heat  pipe  work  as  parallel  fluid  circuits,  with  the 
total  resistance  being  less  when  both  sides  are  heated  versus  applying  heat  to  only  one  side 
for  the  same  heat  flux. 

The  boiling  limit  in  a  flat  heat  pipe  is  similar  to  that  of  a  cylindrical  heat  pipe.  Both 
involve  pressure  balances  at  the  onset  of  nucleation  (pressure  difference  from  bubble  to 
vapor  space  through  liquid  in  the  wick).  However,  since  the  vapor  area  and  heat  input  cross 
section  are  all  rectangular,  the  boiling  limit  equation  for  cylindrical  heat  pipes  (equation  3.1) 
changes  to  equation  4.5  for  flat  heat  pipes. 
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In  equation  4.5,  is  the  vessel  wall  thickness  (m),  is  the  vessel  wall  thermal  conductivity 
(W/m  K),  and  t^  is  the  wick  thickness  (m). 

The  boiling  limit  can  be  expected  to  be  a  very  large  value.  Flat  heat  pipes  are 
constructed  of  high  strength  metals  with  polished  surfaces.  These  metals  feature  very  small 
surface  defect  sizes  and  therefore  minute  nucleation  site  sizes  for  the  onset  of  boiling.  This 
drives  the  maximum  heat  input  for  the  boiling  limit  up  to  a  large  value.  This  limit  will  be 
much  higher  than  the  capillary  limit,  which  will  be  retained  as  the  maximum  heating  rate  for 
the  project  in  consideration.  However,  neither  of  these  limits  will  ever  actually  be  reached, 
as  the  internal  pressure  (due  to  the  temperatures  reached  at  these  heating  limits)  could  cause 
structural  failure  of  the  vessel.  Appendix  2  shows  the  boiling  limit  derivation. 

With  water  as  the  working  fluid,  the  entrainment  and  viscous  limits  are  not  an 
operational  limit  consideration,  as  the  physical  properties  of  water  drive  the  entrainment  and 
viscous  limits  to  very  large  values  within  this  project’s  operating  parameters  (20°C  to  100°C 
condenser  temperature).  Furthermore,  the  entrainment  limit  is  usually  not  a  factor  in  well- 
wicked  heat  pipes,  since  the  wick  retains  the  liquid  to  prevent  entrainment  by  the  vapor. 
Being  difficult  to  calculate,  the  entrainment  limit  is  found  by  balancing  the  shear  force  on 
the  liquid  due  to  the  vapor  flow  with  the  surface  tension  of  the  liquid.  For  most  applications, 
the  entrainment  limit  is  calculated  by  assuming  the  Weber  number  equal  to  one,  where  the 
Weber  number  is  a  ratio  of  vapor  inertial  forces  to  liquid  surface  tension  forces.  One  must 
perform  this  balance  at  the  given  operating  point  of  the  heat  pipe  to  determine  if  entrainment 
is  a  problem  (Peterson,  1994). 

In  the  sonic  limit,  the  vapor  cannot  travel  faster  than  the  local  speed  of  sound  in  a 
constant  area  duct.  The  heat  input  that  can  be  delivered  is  therefore  limited  by  the  maximum 
mass  flow  rate  of  the  vapor.  To  see  the  relatively  simple  derivation,  see  Appendix  3.  The 
sonic  limit  is  found  by  equation  4.6,  where  Tio^j,  is  the  local  temperature  (K)  and  R  is  the  gas 
constant  for  the  vapor  [kJ/(kg  K)]: 


Qs  = 


(4.6) 
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4.2.  FLAT  HEAT  PIPE  STRUCTURAL  CONSIDERATIONS 

A  flat  heat  pipe  has  structural  problems  that  are  not  associated  with  cylindrical  heat 
pipes.  Cylindrical  heat  pipes  are  considered  natural  pressure  vessels  due  to  their  shape,  and 
can  withstand  large  pressure  differences  as  well  as  the  resulting  compressive  or  tensile  forces 
on  their  walls.  Since  heat  pipes  can  exist  at  pressures  lower  than  atmospheric  pressure  when 
not  in  use,  and  can  exceed  atmospheric  pressure  at  elevated,  in-uSe  temperatures,  the 
cylindrical  design  has  clear  advantages.  In  this  case,  one  of  the  flat  heat  pipe’s  advantages-a 
large  surface  area  to  volume  ratio-is  detrimental  to  its  structural  integrity.  The  heat  pipe’s 
surface  area  is  subjected  to  the  pressure  difference  created  between  the  heat  pipe  internal 
pressure  and  the  external  environmental  pressure.  This  may  result  in  a  large  force  acting  on 
the  vessel,  which  can  cause  material  failure.  This  was  clearly  a  design  consideration  that 
needed  to  be  solved. 

The  design  operating  point  for  the  flat  heat  pipes  in  this  project  was  one  atmosphere 
of  differential  pressure,  with  design  case  for  loading  taken  at  a  near  perfect  vacuum.  In  this 
scenario,  the  101.325  kPa  (14.67  psia)  of  atmospheric  pressure  would  be  acting  to  deform 
the  heat  pipe  with  nothing  more  than  structural  supports  to  counteract  the  pressure.  Since 
one  of  the  project  goals  was  to  design  structural  supports  that  did  not  inhibit  two-dimensional 
flow,  any  internal  beams  or  walls  were  ruled  out  as  internal  supports.  This  left  only  pin 
supports  as  an  option.  Not  knowing  how  the  pins  were  to  be  fastened  to  the  vessel  wall  in 
the  final  design,  the  pins  were  modelled  as  both  welded  (rigidly  attached)  and  not  welded  to 
the  vessel  wall.  The  vessel  wall  was  then  modelled  as  a  beam,  with  either  the  side  vessel 
wall  and  a  pin  or  two  pins  acting  as  supports.  This  gave  the  maximum  expected  deflection 
for  the  flat  heat  pipe  of  a  given  thickness,  material,  and  spacer  distance.  A  spreadsheet  was 
developed  to  analyze  deflections  for  many  different  materials,  material  thicknesses,  and  pin 
spacing.  To  view  the  derivations  and  a  sample  spreadsheet,  see  Appendix  4. 

Beam  deflection  was  not  the  only  structural  analysis  conducted.  To  determine 
whether  or  not  deflection  was  the  limiting  factor,  a  column  buckling  analysis  for  pins  of 
varying  lengths,  diameters,  materials,  and  types  of  attachments  also  had  to  be  conducted. 
However,  this  calculation  is  only  necessary  to  reinforce  the  value  obtained  from  the 
deflection  problem,  and  was  performed  solely  on  the  case  for  which  a  solution  was  obtained 
via  the  beam  deflection  problem.  The  equations  necessary  for  this  calculation  are  readily 
found  in  the  Strength  of  Materials  textbook  by  Ferdinand  Singer  (1962). 
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5.0  USNA  FLAT  HEAT  PIPE 

The  corporate  sponsor  funding  this  project  instilled  certain  guidelines  for  the  flat  heat 
pipe  being  constructed.  They  determined  the  dimensions,  which  are  1 .22  m  (48  inches)  long 
by  30.48  cm  (12  inches)  wide  by  1.27  cm  (V2  inch)  thick.  They  also  decided  that  the  vessel 
material,  structural  supports,  and  wicking  material  be  made  of  Monel  R-400.  Due  to  facility 
limitations,  safety  issues,  and  difficulties  in  working  with  liquid  metal  heat  pipes,  the  USNA 
heat  pipe  was  established  as  a  low  temperature  heat  pipe.  Therefore,  distilled  water  was 
chosen  to  be  the  working  fluid,  as  it  has  the  highest  figure  of  merit  of  any  of  the  candidate 
working  fluids  at  the  operating  temperature  (100°C).  Since  effective  thermal  conductivity 
is  not  as  important  an  issue  to  this  design  as  maximizing  capillary  pumping  performance,  it 
was  decided  to  construct  the  USNA  flat  heat  pipe  with  a  composite  wick  consisting  of  two 
Monel  400  screens  of  different  porosity.  The  screen  at  the  liquid- vapor  interface  will  be  120 
X  120  square  mesh  (N=120)  to  provide  a  large  capillary  pumping  head,  and  the  screen  in 
contact  with  the  vessel  wall  will  be  40  x  40  (N=40)  square  mesh  to  provide  a  path  of  low 
flow  impedance.  In  order  to  counteract  the  crushing  force  of  the  atmosphere  when  the  heat 
pipe  is  not  at  operating  temperatures.  Monel  pins  were  used  for  internal  column  supports. 
As  discussed  earlier,  pins  offer  structural  support  without  biasing  the  flow  towards  any 
direction. 

Initially,  the  desired  wall  thickness  and  pin  spacing  diameter  was  not  known.  To 
solve  these  problems,  two  different  spreadsheets  were  used.  The  first  modelled  the  area  in 
between  the  first  pin  from  the  wall.  The  second  the  spacing  between  the  pins.  The 
derivations  of  these  modelling  conditions  and  the  equations  that  were  solved  for  in  terms  of 
pin  spacing  and  material  thickness  are  included  in  Appendix  4  and  are  also  discussed  in 
section  5.2.  Using  a  targeted  maximum  deformation  of  0.15875  cm  (1/16  of  an  inch),  it  was 
found  that  0.127cm  (0.050  inch)  thick  Monel  sheets  would  be  sufficient  with  three  inch  pin 
spacing.  The  three  inch  spacing  is  also  convenient  because  three  is  a  multiple  of  twelve  and 
will  allow  uniform  pin  spacing  throughout  the  entire  pipe. 

Once  the  pin  spacing  was  determined,  it  was  necessary  to  solve  for  pin  diameters. 
Initially,  the  pins  were  only  going  to  be  welded  on  one  side  to  the  vessel  and  wick.  Thus 
they  were  modelled  as  columns  rigidly  attached  on  one  end  and  free  on  the  other.  A  column 
buckling  analysis  was  performed  for  the  maximum  loading  conditions  as  well  as  compressive 
failure.  A  factor  of  safety  was  built  into  the  analysis  since  it  was  assumed  that  the  pins 
existed  in  a  perfect  vacuum  (when  in  reality  there  will  only  exist  the  partial  pressure  of  water 
for  the  given  temperature)  and  support  the  entire  crushing  force  of  the  atmosphere.  It  was 
also  assumed  that  all  of  the  pins  adjacent  to  the  pin  being  analyzed  would  fail,  and  that  it 
would  bear  the  weight  of  the  entire  grid  surrounding  itself  to  the  adjacent  pins.  It  was 
discovered  that  compressive  failure  would  be  the  limiting  factor  for  the  pins,  and  their  proper 
diameter  was  determined  from  this  analysis.  The  minimum  pin  diameter  using  this  method 
was  found  to  be  approximately  0.3556  cm  (0.14  inch).  However,  it  was  not  possible  to 
locate  a  vendor  that  could  supply  rigid  Monel  pins  with  diameters  this  small.  The  smallest 
diameter  readily  available  was  0.635cm  (0.25  inch)  and  this  is  what  was  used  to  create  the 
USNA  flat  heat  pipe.  This  also  provides  a  factor  of  safety  for  failure  due  to  buckling. 
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5.1  USNA  FLAT  HEAT  PIPE  HEATING  LIMITS 

Using  equations  4.1,  4.3,  and  4.4  explained  in  section  5.2  and  derived  in  the 
appendices,  the  capillary  limits  are  as  follows: 

=  426974 

for  the  case  with  only  one  side  heated  and 

ri  =  853914 

le  m 


for  the  case  with  both  sides  heated. 

For  the  case  with  only  one  side  heated,  with  the  given  heater  area,  this  results  in  a 
max  heat  input  of  3967W,  and  7933W  for  the  case  with  both  sides  heated. 

The  boiling  limit,  defined  by  equation  4.5  and  the  derivations  in  Appendix  2,  is 
calculated  to  have  a  maximum  heat  input  of 9527W  assuming  a  typical  nucleation  site  radius 
of  2.5  X  lO"'’  m 


The  sonic  limit  was  calculated  as  325,042W  using  HTPIPE,  a  one-dimensional 
steady-state  heat  limit  prediction  program  for  heat  pipes.  HTPIPE  was  not  accurate  for  the 
other  heating  limits  due  to  approximations  that  one  must  make  to  input  parameters  into  the 
program  which  can  not  accurately  account  for  the  composite  wick  construction. 

There  are  no  standard  methods  of  construction  for  flat  heat  pipes,  as  they  are  a 
relatively  new  and  unique  device.  To  help  remedy  this,  a  subscale  model  was  designed  and 
constructed  to  help  facilitate  the  fabrication  process  with  the  fullscale  model.  Although  not 
necessary,  the  subscale  model  will  have  internal  pin  supports.  These  exist  solely  to  design 
the  final  construction  procedure  for  attaching  these  pins  to  the  vessel  wall  in  a  manner  where 
they  will  provide  the  intended  structural  support.  This  subscale  model  also  has  a  fitting  on 
each  end  to  perform  charging  experiments  and  perfect  the  evacuation  and  charging 
procedure.  The  subscale  model  will  also  be  used  to  calibrate  laboratory  equipment  while 
waiting  for  the  construction  of  the  final  design.  The  subscale  model  is  a  fully  functional  flat 
heat  pipe  and  will  provide  useful  data  to  contribute  to  the  overall  mission  of  the  project. 
Figures  8-11  depict  the  subscale  and  fullscale  heat  pipe  designs. 


Note:  Figures  ore  not  drown  to  scale 
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Figure  9.  Subscale  Heat  Pipe  Interior 
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Note:  Figures  are  not  drawn  to  scale 
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6.0  SUBSCALE  FLAT  HEAT  PIPE 

After  the  necessary  supplies  were  obtained,  the  fabrication  and  assembly  phase  could 
be  carried  out.  The  subscale  heat  pipe  was  constructed  first.  Its  purpose  was  to  provide 
practical  experience  of  flat  heat  pipe  construction  as  well  as  be  used  to  help  calibrate 
instruments  and  set  up  the  laboratory  spaces.  In  addition,  it  aided  in  proving  the  practicality 
of  the  design  and  revealing  and  structural  faults. 

6. 1  SUBSCALE  FLAT  HEAT  PIPE  CONSTRUCTION 

The  basic  approach  to  construction  was  to  assemble  the  heat  pipe  in  two  separate 
components.  This  involved  drilling  all  the  holes  necessary  for  the  pin  supports,  bending  the 
sides  of  the  two  halves  to  the  proper  dimensions,  and  attaching  the  two  layers  of  screen.  The 
end  pieces  and  fittings  were  also  fabricated  at  this  time  and  the  pins  were  cut  to  the  proper 
lengths  and  deburred.  Then,  the  two  sides  were  welded  together,  the  ends  welded  on,  and 
the  fitting  welded  in.  After  this,  the  pins  were  inserted  and  welded  to  both  sides  of  the  heat 
pipe  from  the  outside. 

The  subscale  heat  pipe  was  constructed  mostly  to  the  above  parameters.  The  monel 
sheets  were  bent  using  twelve  inch  long  copper  bars  with  edges  of  fine  radii.  The  screen  was 
then  cut  to  the  same  size  using  metal  shears  and  taped  to  the  vessel  material.  They  were  then 
tack-welded  into  place,  centered  halfway  in  between  pin  locations.  Next,  the  0.635  cm  (1/4 
inch)  holes  for  the  three  spacers  were  carefully  dimensioned  and  drilled  in  both  halves  through 
the  vessel  material  and  the  screen.  This  created  some  problems  with  the  screen  integrity 
around  the  hole  as  well  inducing  stress  around  the  tack  weld  spots.  However,  these  were 
very  minor  and  not  enough  to  substantiate  reattaching  different  screen  material  as  the  screen 
still  maintained  good  pore  integrity  and  contact  with  the  vessel  wall.  This  meant  that  the  wick 
would  still  produce  sufficient  capillary  pumping  action  and  the  effective  thermal  conductivity 
of  the  screen  would  remain  as  high  as  possible.  The  end  fittings  for  the  high  vacuum  ball 
valves  were  then  machined  and  the  endpieces  and  the  pin  spacers  were  then  cut.  However, 
they  required  extensive  deburring  and  grinding  so  that  they  could  fit  into  the  0.635  cm  (1/4 
inch)  diameter  spacer  hole. 

After  the  subscale  flat  heat  pipe’s  constituent  elements  were  readied,  the  two  halves 
were  welded  together.  The  welds  were  staged  from  side  to  side  so  that  the  warping  due  to 
thermal  stress  would  be  minimized.  The  ends  were  then  welded  on,  and  the  fittings  were 
attached  in  the  proper  place.  Following  this,  the  pins  were  Tungsten  Inert  Gas  (TIG)  welded 
to  both  sides  of  the  vessel  wall  from  the  outside. 

6.2  SUBSCALE  HEAT  PIPE  TESTING 

Once  completed,  the  subscale  was  pressure  tested  for  leaks.  The  pins  presented  a 
problem,  as  did  the  end  fittings.  To  make  an  airtight  seal,  the  end  fitting  design  had  to  be 
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alterred.  The  tubes  that  extended  from  the  heat  pipe  were  cut,  and  the  fittings  were  welded 
directly  to  the  ends  of  the  heat  pipe.  The  pins  were  welded  to  both  sides  of  the  vessel,  as  it 
would  have  been  very  difficult  to  attach  the  pins  to  only  one  side  of  the  heat  pipe.  This 
reduced  the  calculated  deflection  substantially,  and  would  later  do  the  same  for  the  final 
model. 


The  pins  required  re-welding  several  times  in  order  to  create  an  airtight  seal  at  the 
desired  1  atmosphere  of  gage  pressure.  This  problem  of  rewelding  was  found  to  be  inherent 
to  the  design  by  both  welder  and  design  team.  The  pins  were  necessary  for  structural  integrity 
and  had  to  be  welded  from  both  sides  to  ensure  structural  support.  Reheating  of  a  completed 
weld  by  welding  on  the  opposing  side  was  determined  not  to  be  a  pertinent  issue  since  the 
type  of  welding  utilized  featured  the  use  of  localized  heating  for  a  short  period  of  time.  This 
would  not  be  sufficient  to  heat  the  opposing  weld  spot  to  a  temperature  at  which  it  could  start 
to  loosen  (Monel  welding  wire  or  sheets  would  liquefy). 

The  subscale  heat  pipe  was  not  intended  to  be  used  for  extensive  experimental  data- 
taking,  but  to  be  used  to  verify  that  the  basic  heat  pipe  design  worked,  to  find  design  flaws, 
and  to  setup  laboratory  equipment.  Therefore,  its  testing  regime  was  not  very  expansive.  Its 
testing  regime  and  objectives  were  as  follows: 

1 .  Design  a  reliable,  easily  repeatable  method  of  leak-testing  and  make  the  subscale  heat  pipe 
leak-free. 

2.  Evacuate  the  heat  pipe.  Attach  a  pressure  gauge  to  the  evacuation  apparatus  and  measure 
the  level  of  evacuation  achieved. 

3 .  After  successful  evacuation,  design  and  construct  a  charging  apparatus  for  accurately 
charging  the  heat  pipe  with  the  correct  amount  of  working  fluid  and  sealing  it  off  after 
charging.  This  method  must  be  easily  repeatable. 

4.  Apply  a  heat  source  to  one  end  of  the  heat  pipe  and  check  for  response.  The  condenser 
end  will  achieve  a  temperature  nearly  the  same  temperature  as  the  end  with  the  applied  heat 
source.  It  is  advisable  not  to  use  anything  more  than  hot  water  for  this  test  as  it  is  conducted 
by  hand. 

5.  Submerge  one  end  of  the  heat  pipe  in  a  hot  water  bath.  Monitor  the  water  bath 
temperature  and  the  heat  pipe  by  quantifiable  means.  Also  monitor  the  heat  pipe  performance 
with  the  infrared  camera  to  determine  if  it  is  working  properly,  checking  for  end  effects  that 
might  indicate  a  leak  or  improper  charging  and  temperature  variations  that  signal  interior 
defects  or  improper  charging. 


6.3  SUBSCALE  HEAT  PIPE  TESTING  RESULTS 

The  leak-testing  and  charging  consisted  of  the  fabrication  of  a  charging  apparatus. 
A  0.635  cm  (1/4  inch)  nipple  was  fitted  to  the  end  of  it  (with  necessary  reducer)  and  mated 
to  the  compressor.  The  heat  pipe  was  then  pressurized  for  leak  testing.  Valve  2  in  the 


34 


charging  apparatus  (see  Figure  12)  could  be  throttled  so  that  excessive  pressure  did  not  build 
up  in  the  heat  pipe.  A  pressure  gauge  was  fitted  temporarily  to  the  charging  apparatus  near 
the  heat  pipe  to  calibrate  the  throttling  and  was  removed  at  the  time  of  this  picture.  To  leak 
test,  soapy  water  was  applied  to  the  pressurized  heat  pipe  in  order  to  detect  the  leaks.  Leaks 
would  form  bubbles  in  the  soapy  water.  This  assembly  could  be  mounted  and  removed  from 
the  heat  pipe  with  a  few  hand  tools. 


Figure  12.  Charging  Apparatus 


Charging  the  heat  pipe  with  working  fluid  was  performed  fairly  simply  as  well.  To 
charge  a  heat  pipe  the  wick  must  be  completely  saturated.  First,  the  heat  pipe  and  charging 
assembly  were  mated  to  the  vacuum  pump  and  evacuated  with  valve  2  closed.  In  the  fluid 
charging  column  was  placed  the  correct  amount  of  water  to  charge  the  heat  pipe.  These 
amounts  were  measured  and  marked  (observe  gradations  on  fluid  charging  column),  taking 
into  account  the  volume  that  would  occupy  the  fittings  as  well  as  the  heat  pipe.  Once 
evacuation  was  complete,  valve  1  and  the  valve  attached  to  the  heat  pipe  (not  pictured)  were 
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closed,  and  then  the  vacuum  pump  was  shut  off.  Valve  2  was  then  opened,  and  the  vacuum 
inside  the  fittings  drew  the  water  in  to  fully  fill  the  fittings.  The  heat  pipe  valve  was  then 
opened  slightly  to  bleed  in  the  necessary  charge  (as  marked  on  the  column).  All  valves  were 
then  closed,  and  the  charging  apparatus  removed  if  necessary,  taking  care  not  to  disturb  the 
valve  connected  to  the  heat  pipe  (as  it  holds  the  heat  pipe  to  the  proper  internal  conditions). 
It  is  important  to  note  that  all  valves  used  were  rated  for  special  high-vacuum  service.  To 
verify  the  time  required  for  charging  the  heat  pipe,  the  mass  flow  rate  at  6.9  kPa  (Ipsia)  was 
calculated  for  water  vapor  through  the  0.635  cm  (1/4  inch)  hole  in  the  end  fitting  through 
which  the  pipe  was  evacuated.  The  vacuum  pump’s  flow  rate  was  found  to  be  the  limiting 
factor  in  evacuation,  not  choked  flow  conditions.  This  condition  would  remain  true  for  the 
fullscale  heat  pipe  as  it  was  to  have  the  same  size  end  fitting. 

The  finished,  testable  model  of  the  subscale  flat  heat  pipe  appears  in  Figure  13.  It  is 
shown  with  a  valve  still  attached  to  it  to  show  how  they  are  mounted.  Note  the  absence  of 
the  tube  in  the  end  fitting.  The  heat  pipe  is  painted  flat  black  to  raise  its  surface  emissivity 
so  it  can  be  viewed  best  by  an  infrared  camera. 


Figure  13.  Subscale  Heat  Pipe 


After  all  leaks  were  eliminated,  and  the  heat  pipe  was  successfully  charged,  it  was 
subjected  to  a  hot  water  bath.  Within  seconds,  the  other  side  of  the  heat  pipe  began  to  warm 
up.  It  worked  at  least  as  a  crude  heat  pipe.  After  this,  it  was  submerged  in  a  pot  of  boiling 
water  and  monitored  with  the  infrared  camera.  This  facilitated  a  qualitative  check  of  the  heat 
pipe  performance  as  well  as  calibrating  and  learning  how  to  use  the  infrared  camera.  The 
subscale  passed  this  test  ftawiessly,  showing  an  almost  isothermal  condenser  section  with 
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marginal  end  effects  (slightly  cooler  temperatures)  due  to  the  fin  effect  of  the  end  fittings  and 
the  innate  properties  of  a  rectangular  object  (corners  and  ends  are  cooler  than  the  middle 
when  heated).  See  Figure  14  for  a  photo  of  the  infrared  camera. 


Figure  14.  Infrared  Camera 


6,4  FULLSCALE  HEAT  PIPE  IMPLICATIONS 

After  testing  the  subscale  model,  it  was  noted  that  several  features  initially  designed 
into  the  fullscale  model  were  no  longer  necessary. 

Only  one  fitting  was  found  necessary  to  perform  complete  leak  testing,  evacuation, 
and  charging.  Therefore,  the  fullscale  has  only  one  end  fitting.  The  tube  in  the  end  fitting  was 
weak  and  unnecessary,  thus  it  was  eliminated.  Not  only  is  it  easier  from  a  fabrication 
standpoint,  as  it  aids  alignment  from  side  to  side,  but  it  also  provides  greater  structural 
support.  Thus,  the  pins  were  welded  on  both  sides  to  the  vessel  (as  was  done  in  the  subscale 
model).  The  shop  received  new  equipment  between  construction  of  the  two  heat  pipes,  and 
this  affected  the  constaiction  process.  See  Section  7  for  detail  of  the  fullscale  constniction 
process. 
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7,0  FULLSCALE  HEAT  PIPE  CONSTRUCTION 

The  flillscale  model  was  to  be  constaicted  differently  than  the  subscale  model.  This 
was  due  to  Technical  Support  Department's  acquisition  of  a  new  CNC  Punch-Press  machine 
as  well  as  experimentation  with  the  subscale  and  a  visit  from  the  corporate  sponsor. 

7.1  FULLSCALE  HEAT  PIPE  CONSTRUCTION  DETAILS 

After  the  screens  and  Monel  sheets  were  cut  to  the  proper  dimensions,  their  size  and 
assigned  hole  placement  were  loaded  into  the  CNC  Punch-Press  machine's  user  interface.  The 
Monel  sheets  for  the  vessel  wall  were  punched  first,  making  holes  in  the  locations  where  pins 
were  to  be  welded.  Next,  to  avoid  any  tearing,  stretching,  or  any  other  kind  of  deformation 
in  the  screen,  the  screens  were  laid  on  top  of  each  to  other,  taped  to  the  side  of  thin  steel 
sheets  (30mil),  and  sandwiched  in  between  the  sheets.  They  were  then  loaded  onto  the  CNC 
Punch-Press  machine  and  were  punched  in  the  appropriate  places,  so  that  the  steel  sheets 
absorb  the  initial  impact  and  prevent  the  movement  and  subsequent  deformation  of  the  screen. 
By  using  the  computerized  machine,  it  was  possible  to  place  all  of  the  necessary  holes  in 
exactly  the  correct  spots  on  both  the  sheets  and  the  screen  to  ensure  that  the  sides  would  line 
up  perfectly. 

Long  copper  bars  with  fine  radius  edges  (same  radius  as  used  in  construction  of 
subscale)  were  obtained  to  facilitate  bending  the  flillscale  heat  pipe  to  the  required 
dimensions.  The  two  layers  of  screen  were  then  layed  on  top  of  the  vessel,  fastened  with  tape 
on  the  ends,  and  fastened  down  to  the  vessel  by  threading  small  screws  through  the  pin  spacer 
holes  and  attaching  nuts  on  either  sides.  The  screen  was  then  tack- welded  to  the  vessel  wall 
in  regular  intervals  between  the  pin  spacer  locations.  Keeping  the  screen  flush  and  in  contact 
with  the  vessel  presented  a  challenge,  and  a  small  portion  of  the  screen  on  one  side  of  the  heat 
pipe  ended  up  being  not  in  contact  with  the  vessel,  but  detached  from  the  vessel  in  a  small 
section  of  the  heat  pipe.  This  spot  was  noted  in  case  it  caused  any  performance  problems. 
Figure  1 5  shows  a  close-up  photo  of  the  tack  welds  and  the  holes  drilled  for  pins  in  a  Monel 
sheet,  and  Figure  16  shows  the  two  completed  halves  lying  on  top  of  each  other.  The  end 
pieces  and  fitting  were  machined.  Figure  17  shows  a  completed  half  with  visible  screen  and 
bent  side.  The  pins  were  then  cut  to  the  proper  length,  milled,  and  deburred  to  fit  into  the 
necessary  space.  Figure  18  shows  the  pins  standing  on  the  vessel  wall  before  this  process. 

With  all  of  the  necessary  components  ready,  the  assembly  began.  First,  the  sides  of 
the  two  separate  halves  were  welded  together,  carefully  sequencing  the  welds  and  using  a 
heat  trap  to  minimize  deformation.  After  this,  the  ends  were  welded  on,  and  the  singular  end 
fitting  was  welded  snug  to  the  end  to  resemble  the  final  appearance  of  the  subscale  heat  pipe. 
Following  this,  the  pins  were  placed  in  their  proper  spots  and  TIG  welded  from  either  side 
of  the  heat  pipe.  The  experimental  setup  then  began. 
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8.0  EXPERIMENTAL  SETUP  AND  PROCEDURE 

The  setup  for  experimentation  as  well  as  the  procedure  for  experimentation  were 
complex  parts  of  the  project  and  had  great  importance  to  the  success  of  the  project.  It  was 
very  important  to  verify  that  the  flat  heat  pipe  does  indeed  perform  as  a  heat  pipe  and  to  take 
data  on  its  one-dimensional  performance  before  proceeding  on  to  any  other  tests.  These  one¬ 
dimensional  tests  were  also  the  easiest  way  to  understand  and  interpret  the  internal  dynamics 
of  the  heat  pipe  as  they  influenced  the  temperature  profile.  Next,  as  the  impetus  behind  the 
project  was  to  examine  two-dimensional  flat  heat  pipe  behavior,  asymmetric  heating 
conditions  were  imposed  and  the  experiments  run.  A  heat  pipe  is,  by  definition,  a  hermetically 
sealed  vessel  that  transfers  heat  via  phase  change  processes  and  has  capillary  pumping  action 
to  circulate  the  working  fluid.  It  was  necessary  to  test  this  capillary  action.  As  mentioned 
earlier,  the  capillary  limit  is  also  the  desired  operating  limit  for  most  heat  pipes.  Finding  the 
capillary  limit  of  the  heat  pipe  is  therefore  very  important  to  not  only  understanding  heat  pipe 
performance  but  also  to  maximizing  heat  pipe  performance.  Experiments  were  conducted 
at  differing  angles  of  gravity  assist  for  the  liquid  return  path  to  test  the  capillary  pumping 
action  of  the  screen  and  gravity’s  effect  on  the  heat  pipe’s  performance. 

8.1  EXPERIMENTAL  SETUP 

The  asymmetric  heating  case  was  a  unique  case.  The  heaters  had  to  possess  the 
flexibihty  to  be  reconfigured  and  turned  on  and  off  independently  to  be  able  to  satisfy  both 
the  symmetric  and  asymmetric  heating  cases.  To  meet  this  need,  16  high  temperature  500 
Watt  strip  heaters  1.905  cm  (3/4  inch)  inch  wide  with  a  heated  section  of  30  cm  (11  13/16 
inches)  were  purchased.  They  were  wired  in  blocks  of  four  to  variable  transformers,  and 
could  therefore  be  moved  in  blocks  of  four  and  have  their  heating  rate  adjusted  independently 
in  four  blocks  of  four  heaters  each.  The  variable  transformers  were  also  wired  to  Wattmeters, 
which  were  wired  into  junction  boxes  (to  be  detailed  later).  The  electric  signals  fi'om  the 
Wattmeters  were  recorded  by  the  computer  to  which  these  junction  boxes  fed. 

Thermal  joint  compound  was  spread  between  the  heaters  and  the  heat  pipe  to  lower 
contact  resistance  and  ensure  that  most  of  the  generated  heat  went  into  the  heat  pipe.  The 
heat  pipe  was  then  insulated  with  a  ceramic  tile  over  the  heater  sections  with  aluminum  foil 
insulation  wrapped  around  the  entire  assembly.  To  view  the  heaters  being  mounted  with 
thermal  joint  compound  in  the  symmetric  case  (all  heaters  on  one  side,  covering  the  entire 
width)  see  Figure  19.  Figure  19  does  not  show  all  of  the  heaters  mounted  to  the  heat  pipe 
for  the  symmetric  cases.  Instead,  Figure  1 9  shows  ten  of  the  heaters  mounted  and  the  thermal 
joint  compound  that  was  part  of  the  heater  section  fabrication.  Figure  20  shows  the 
completed  asymmetric  mounting  (half  of  the  heaters  on  each  side,  covering  only  half  of  one 
end  with  the  heaters  covered  with  a  ceramic  tile  and  aluminum  foil  insulation). 
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The  same  evacuation  and  charging  assembly  that  was  developed  for  the  subscale  heat 
pipe  was  used  for  the  fullscale  heat  pipe.  A  photo  of  this  assembly  is  shown  in  section  6.  To 
record  saturation  pressures  inside  of  the  heat  pipe,  a  pressure  transducer  was  mounted  to  the 
end  of  the  heat  pipe  in  between  the  heat  pipe  valve  and  the  heat  pipe  itself  Figure  21  shows 
this  end  fitting  with  the  transducer  mounted.  This  transducer  was  used  to  record  most  of  the 
data.  However,  this  transducer  malfunctioned,  and  was  replaced  by  an  analog  pressure  gauge 
as  seen  in  Figure  22. 


Figure  2 1 .  Heat  Pipe  End  Fittings 


Figure  22,  Analog  Pressure  Gauge 


Again,  it  is  important  to  note  that  the  valve  used  was  specially  rated  for  high  vacuum 
service.  A  regular  valve  would  not  be  capable  of  holding  the  high  vacuum  states  that  occur 
in  evacuation  and  when  the  heat  pipe  is  properly  charged  but  not  operating. 

Since  testing  at  various  angles  to  the  horizontal  is  an  important  part  of  the  test  criteria, 
it  was  necessary  to  design  a  test  stand  such  that  the  heat  pipe  was  capable  of  rotating  through 
the  required  angles  while  remaining  free  to  natural  convection.  For  this,  a  wooden  stand  was 
designed  that  consisted  of  two  legs  fixed  to  a  flat  base.  See  Figure  23  to  view  the  stand  with 
heat  pipe  mounted.  Attached  to  the  stand  were  two  clamps  to  hold  the  heat  pipe  in  place 
while  being  able  to  rotate  it.  These  were  made  from  two  C-Clamps  welded  together  by  a  steel 
bar.  From  this  bar  protruded  a  threaded  screw,  which  was  fixed  to  the  bar.  This  screw  went 
through  a  wooden  bearing  block  and  the  test  stand  legs,  and  was  attached  on  the  other  side 
of  the  leg  by  a  wing  nut  and  a  washer.  The  wing  nut  could  be  loosened  and  tightened  to 
move  the  clamps  through  the  various  angles  about  the  screw.  This  assembly  rotated  the  heat 
pipe  through  the  required  test  angles.  Figure  24  shows  the  clamps  mounted  to  the  heat  pipe 
with  PVC  feet  attached.  These  PVC  feet  served  not  only  to  spread  the  load,  but  also  to 
insulate  the  clamps  from  the  heat  pipe  so  that  the  clamps  could  not  become  a  heat  sink. 

Temperature  data  were  recorded  in  two  ways.  The  first  was  via  the  infrared  camera 
videography,  pictured  in  section  6.  The  other  means  was  by  an  array  of  thermocouples.  The 
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placement  of  the  thermocouples  on  the  heat  pipe  are  shown  in  Figure  25.  Figure  26  shows 
how  the  thermocouples  actually  looked  when  mounted  on  the  heat  pipe.  The  numbers  of  the 
thermocouple  locations  correspond  to  their  assigned  numbers  in  the  data  collection  software. 
The  thermocouples  were  wired  into  the  same  data  acquisition  junction  boxes  used  by  the 
Wattmeters.  There  were  8  junction  boxes  with  8  channels  each  for  a  total  of  64  channels  of 
data  sampling.  However,  only  43  were  used  (39  thermocouples  plus  four  Wattmeters).  To 
see  how  the  system  is  configured,  see  Figure  27. 


Figure  23.  Heat  Pipe  Test  Stand 
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Figure  24.  Heat  Pipe  Mounted  to  Clamp 
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8.2  EXPERIMENTAL  PROCEDURE 

The  heat  pipe  went  through  tests  at  different  orientations.  Figure  28  shows  the 
nomenclature  for  the  different  orientations. 


Heater  End 


Figure  28.  Testing  Orientations 


The  testing  criteria  for  the  symmetric  heating  case  (all  heaters  mounted  on  one  side, 
covering  entire  width  of  heat  pipe)  was  as  follows; 

1.  With  100%  fluid  charge,  test  orientations  2,3,4,5,6,  and  7  heating  slowly  to  at  least  400 
Watts  (more  if  conditions  allow),  allowing  heat  pipe  to  reach  equilibrium  temperature  for  each 
heat  input. 

2.  With  125%  fluid  charge,  test  orientations  3,4,5,  and  7  heating  with  at  least  400W  allowing 
equilibrium  conditions. 
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3.  With  75%  charge,  test  orientations  3,4,  and  5  heating  with  at  least  400W  allowing 
equilibrium  conditions. 

4.  Mounting  the  heaters  in  the  asymmetric  configuration  (see  Figure  20),  test  heat  pipe  with 
100%  fluid  charge  in  orientations  3,4,and  5  heating  slowly  to  as  high  a  heat  input  as  possible, 
allowing  for  equilibrium  between  heating  increases. 

5.  For  all  runs,  data  is  to  be  collected  over  time  by  data  acquisition  software,  which  records 
continuously,  and  the  infrared  camera,  which  records  at  the  user’s  discretion. 

6.  Take  video  and  still  photographs  of  interesting  phenomena  via  the  infrared  camera. 
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9.0  DATA  ANALYSIS 


The  experiments  detailed  in  chapter  8.0  were  completed.  However,  some  calculations 
had  to  be  performed  before  the  data  could  be  properly  analyzed. 


A  fin,  whether  it  is  an  isothermal  heat  pipe  fin  or  a  conventional  fin,  rejects  heat  to  the 
ambient  by  means  of  convection  and  thermal  radiation.  These  heat  rejection  modes  are 
represented  by  equations  9.1  and  9.2,  respectively  (Kreith  and  Brohn,  1997).  In  equation  9.1, 
h  is  the  convection  coefficient  [W/(m  K)],  A  is  the  surface  area  of  the  body,  and  AT  is  the 
temperature  difference  between  the  hot  body  and  its  surroundings.  In  equation  9.2,  a  is  the 
Stefan-Boltzmann  constant,  s  is  surface  emissivity  (percent  radiated  as  compared  to  an  ideal 
emitter),  T,,  is  the  temperature  of  the  radiating  body,  and  T^^^  is  the  ambient  environment 
temperature. 

(9.1) 


Qcon.  =  hA^T 


(9.2) 

a-  =  -  tl) 


For  radiation  heat  transfer,  the  emissivity  is  set  to  .95,  a  value  experimentally 
determined  for  the  heat  pipe  using  the  infrared  camera  and  thermocouple  measurements.  This 
was  done  by  comparison  of  thermocouple  measurement  to  injfrared  measurement.  Emissivity 
is  adjusted  in  the  infrared  camera  data  menu  until  the  temperatures  seen  by  the  infrared 
camera  match  the  temperatures  recorded  by  the  thermocouples. 

Tables  and  equations  used  in  solving  the  combined  convection-radiation  heat  rejection 
problem  are  shown  in  Appendix  5.  Figure  29  shows  the  heat  transfer  rejection  modes  for  a 
given  heat  input  (Watts)  assuming  all  heat  is  transferred  into  the  heat  pipe  (Cerza,  1999). 

From  Figure  29  it  is  seen  that  both  radiation  and  convection  are  important  modes  of 
heat  transfer  for  the  heat  pipe.  To  determine  if  the  heat  pipe  was  properly  functioning  as  a 
heat  pipe,  it  can  be  compared  to  a  conventional  fin.  A  fin  is  an  extended  surface  that  enhances 
heat  transfer  by  increasing  the  heat  transfer  surface  area,  such  as  the  aluminum  fins  on  an 
automobile  radiator.  Thermal  energy  is  input  to  the  fin  at  one  end.  The  energy  then  conducts 
through  the  fin  to  the  colder  end.  As  the  energy  conducts  through  the  fin  it  is  rejected  to  the 
ambient  surroundings  by  convection  and  radiation  (in  this  case).  Since  convection  and 
radiation  are  of  the  same  order  of  magnitude  for  the  heat  input  range  of  these  experiments, 
the  simpler  convection  model  is  used  with  a  convection  heat  transfer  coefficient,  h,  equal  to 
twice  that  of  the  value  for  convection  only,  thus  taking  the  radiation  heat  rejection  mode  into 
account.  Steady-state  fin  temperature  along  the  fin  length  can  then  be  calculated  for  a 
rectangular  fin  by  equation  9.3. 
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~  convection 
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Figure  29,  Heat  Transfer  Rejection  Modes  for  Heat  Input 
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where  m  is  found  with  equation  9.3  as  follows: 


111  = 


f— 1 

V  ^  j 


(9.3) 


(9.4) 


and  P  is  the  perimeter  of  the  fin,  is  the  cross-sectional  area  in  the  x-direction. 

The  previous  fin  analysis  is  for  a  conventional  fin.  The  heat  pipe  fin  temperature  is 
approximately  isothermal  along  the  fin  length.  Equation  9.5  represents  the  steady-state, 
isothermal  heat  pipe  fin  relation: 

(9.5) 


Equation  9,5  is  a  fourth  order  polynomial  equation  for  the  temperature  T.  The 
solution  to  equation  9.5  is  the  positive  real  root  of  equation  9,5  assuming  that  T  is  unique  and 
greater  than  T.^,,,. 
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9,1  HEAT  PIPE  VERIFICATION 

Figure  30  shows  a  plot  of  a  conventional  fin  vs  the  experimental  data  for  the  USNA 
flat  heat  pipe  at  a  400W  steady-state  condition  with  100%  fluid  charge. 
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Figure  30.  Standard  Fin  vs  Ideal  Heat  Pipe  Fin  and  Actual  Heat  Pipe  (400  W 

Steady-State,  Horizontal) 


As  one  can  see  from  the  experimental  data  plot,  the  heat  pipe  was  not  isothermal,  but 
had  about  a  60°C  temperature  drop  from  the  beginning  of  the  heat  pipe  to  the  end.  However, 
a  portion  of  the  heat  pipe  was  nearly  isothermal.  It  can  be  assumed  that  the  heat  pipe  is 
peiforming  as  a  heat  pipe  for  the  first  80  cm  away  from  the  heater  section  and  not  as  a 
con  ventional  fin  since  the  temperatures  of  the  experimental  data  are  close  compared  to  the 
theoretical  heat  pipe  case.  Figure  3 1  shows  an  infrared  camera  photograph  of  the  heat  pipe 
in  the  horizontal  (3)  position  (not  yet  reached  steady-state  with  125%  charge).  This  clearly 
shows  that  although  there  are  some  serious  end  effects  involved,  a  large  region  of  the  heat 
pipe  is  fairly  isothermal.  It  is  believed  that  some  air  leaked  into  the  heat  pipe  and  caused  the 
end  effects  (cooler  temperature  pockets)  by  blocking  a  portion  of  the  condenser  region.  At 
higher  heat  inputs  (such  as  pictured  in  Figure  3 1 ),  the  air  would  be  compressed  into  the  end 
fitting  by  the  water  vapor,  causing  the  cold  region  to  recede,  making  the  heat  pipe  surface 
closer  to  isothermal  The  cold  spot  in  the  middle  of  the  photo  corresponds  to  a  section  of 
epoxy  used  to  seal  a  leak  It  is  an  insulator,  and  is  therefore  colder  than  the  metal  of  the  heat 
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pipe.  The  clamps  that  hold  the  heat  pipe  to  the  test  stand  also  insulate  and  therefore  are 
shown  to  be  cooler  than  the  heat  pipe.  Figure  32  is  a  graph  corresponding  to  the  photo  in 
Figure  3 1,  the  best  isothermal  profile  achieved  by  the  USNA  flat  heat  pipe. 


Figure  31.  IR  Photo  of  Horizontal  Heat  Pipe,  800W,  125%  Charge 
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Figure  32.  Near  Isothermal  Case  Temperature  Vs  Time,  Florizontal,  600W,  100%  Charge 
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9.2  EFFECTS  OF  ORIENTATION  ON  HEAT  PIPES 

The  convection  coefficient,  h,  utilized  in  equation  9.5  change  for  flat  plates  with  their 
heated  side  facing  downwards  in  accordance  with  their  orientation.  Equations  9.6  and  9.7 
can  be  used  to  solve  for  the  convection  coefficient  at  a  given  angle  of  orientation,  0,  from  the 
vertical.  In  equation  9.6,  Nu  is  the  Nusselt  number,  a  dimensionless  heat  transfer  coefficient 
that  is  the  ratio  of  convection  heat  transfer  to  conduction  in  a  fluid  layer  of  thickness  L.  Gq 
and  Pr  are  the  Grashof  and  Prandtl  numbers,  respectively.  The  Grashof  number  is  the  ratio 
of  buoyancy  to  viscous  forces,  and  the  Prandtl  number  is  the  ratio  of  molecular  momentum 
difflisivity  to  thermal  diffusivity  (Kreith  and  Brohn,  1997). 


(9.6) 


h  = 


km 

L 


=  .56(G/}Prcos^)" 


(9.7) 


Note;  equation  9.7  is  only  valid  for  0<9O° 

For  the  90°  case,  equation  9.8  is  to  be  used,  where  Ra,  is  the  Rayleigh  number, 
which  is  the  product  of  the  Grashof  and  Prandtl  numbers  (BCreith  and  Brohn,  1997). 


Nu=15{Ra,)" 


(9.8) 


Although  the  convection  coefficient  would  be  slightly  higher  for  the  heated  side  facing 
upwards  due  to  the  interaction  of  buoyancy  and  momentum  effects  in  natural  convection,  this 
case  will  be  treated  as  if  the  convection  coefficient  is  uniform  for  both  sides  of  the  heat  pipe. 
By  looking  at  equation  9.7,  one  can  see  that  convection  coefficients  decline  as  the  angle 
increases.  Therefore,  at  orientations  approaching  the  horizontal,  to  reject  the  same  amount 
of  heat,  the  surface  temperature  of  the  heat  pipe  must  rise  (equation  9.1).  Using  equation  9.7, 
a  theoretical  plot  of  isothermal  flat  heat  pipe  temperature  versus  angle  orientation  was  made 
and  compared  against  experimental  data  collected  for  the  400W  steady-state  condition  with 
100%  charge  (Figure  33).  The  experimental  data  is  relatively  close  to  the  theoretical  values. 

When  an  orientation  angle  other  than  horizontal  was  introduced,  the  heat  pipe  did  not 
demonstrate  isothermal  behavior  until  the  heat  input  was  increased  to  near  600W.  Figure  34 
shows  the  best  isothermal  profile  reached  for  the  vertical  orientation  (position  5).  Figure  34 
depicts  the  sides  of  the  heat  pipe  as  much  hotter  than  the  middle  portion.  This  condition 
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existed  for  all  inclined  cases.  Figure  34  shows  a  temperature  difference  of  about  27°C  from 
the  hottest  part  to  the  coolest  part  of  the  condenser  section.  The  data  corresponding  to 
Figure  34  was  the  closest  to  a  theoretical  heat  pipe  of  any  of  the  inclined  test  cases.  Figure 
35  shows  a  start-up  condition  (transient  effect)  for  this  orientation  (position  5)  at  a  lower 
power  setting  of  400W.  Note  the  large  temperature  gradient  across  the  heat  pipe. 


Figure  33.  Heat  Pipe  Temperature  vs  Angle,  400W 
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Figure  34.  IR  Photograph,  Vertical  (5),  600W 


Figure  35.  Vertical  Gradient,  400W,  125%  Charge 
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Only  one  case  was  tested  with  the  condenser  end  lower  than  the  evaporator  end.  This 
was  orientation  2.  The  heat  pipe  did  not  function  properly,  acting  as  a  fin.  This  is  because 
the  maximum  vertical  capillary  pumping  height  achievable  by  the  wick  is  around  five  cm. 
Therefore,  no  water  could  be  delivered  to  the  evaporator  region  via  the  wick  (its  height  was 
greater  than  five  cm  above  the  end  of  the  heat  pipe)  and  the  evaporator  section  remained  dry. 
Without  the  circulation  of  the  working  fluid,  the  heat  pipe  does  not  function  as  a  phase 
change  device,  but  acts  as  a  conventional  fin.  Only  the  horizontal  tests  performed  well,  and 
this  will  be  explained  in  section  9.6. 

9.3  EFFECTS  OF  VARYING  CHARGE 

To  compare  the  effects  of  varying  charges,  each  charge  experiment  was  conducted 
at  a  400W  heat  input.  The  experiments  were  allowed  to  reach  a  steady-state  condition  for 
orientations  3,  4,  and  5.  The  horizontal  cases  produced  the  best  data  and  are  presented  for 
comparison.  Figures  36, 37,  and  38  show  the  data  plots  of  the  each  case,  125%  charge,  100% 
charge,  and  75%  charge  in  sequence.  Figure  39  depicts  the  three  charges  plotted  against  each 
other  on  the  same  graph  for  a  uniform  heat  input  over  time. 


time  (min) 


Figure  36.  Temperature  vs  Time.  400W  125%  Charge 
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Figure  39.  Varying  Charges  Temperature  Vs.  Time,  400W,  Horizontal 


Although  the  differences  in  effects  are  not  very  pronounced,  it  is  apparent  that  the 
75%  charge  case  displays  the  quickest  response  time.  This  is  because  there  is  less  water  in 
the  heat  pipe  and  therefore  less  energy  to  be  absorbed,  allowing  the  heat  pipe  to  attain  steady 
state.  However,  notice  that  not  all  of  the  heat  pipe  in  the  75%  charging  case  achieves  the 
same  operating  temperatures.  This  could  be  due  to  a  lack  of  fluid  in  the  wick  required  to 
circulate  the  necessary  heat  to  activate  the  full  heat  pipe  or  the  presence  of  air  in  the  heat  pipe 
(gas  loading). 

The  1 00%  charging  case  has  the  best  traditional  heat  pipe  response.  When  a  location 
encounters  the  condensing  vapor  front,  it  changes  temperature  rapidly  to  the  saturation 
temperature  corresponding  to  the  internal  saturation  pressure.  However,  it  takes  longer  for 
the  1 00%  charging  case  to  reach  isothermal  conditions  along  the  length  of  the  heat  pipe  than 
the  125%  case.  This  is  due  to  the  presence  of  air,  as  once  it  is  compressed  flirther  towards 
the  end  of  the  heat  pipe  the  temperature  at  the  location  it  previously  occupied  increasing 
rapidly.  Also,  the  100%  charging  case  achieves  the  best  near-isothermal  condition,  as  all 
locations  except  the  end  are  within  a  few  degrees  of  each  other  at  the  steady-state  condition. 
The  divergent  line  on  the  top  of  the  plot,  corresponding  to  1 1 .4cm,  does  not  necessarily  make 
the  heat  pipe  non-isothermal.  This  curve  corresponds  to  thermocouples  on  the  back  of  the 
heater  side,  and  are  therefore  more  representative  of  the  heater  temperatures.  The 
thermocouples  were  reading  an  averaged  temperature  of  the  hot  heaters  and  the  inside 
evaporator  section  of  the  heat  pipe.  Hence,  they  would  indicate  higher  temperatures  than  the 
thermocouples  on  the  condenser  end  of  the  heat  pipe.  This  curve  can  therefore  be  ignored. 
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The  125%  charging  case  theoretically  should  take  more  time  for  all  locations  to  reach 
an  isothermal  condition.  This  is  because  there  may  be  some  liquid  flooding  in  the  evaporator 
section,  and  until  it  is  vaporized,  there  can  be  no  menisci  in  the  wick  and  therefore  no 
capillary  pumping  action.  The  125%  charging  case  does  reach  a  near  isothermal  case  when 
it  reaches  steady-state,  though  it  heats  up  faster  than  the  1 00%  charge.  This  is  not  fully 
understood,  though  it  is  likely  that  the  100%  charging  case  was  conducted  with  more  air 
present  in  the  heat  pipe  than  the  125%  charging  case,  which  required  more  time  to  compress 
the  air  towards  the  end  of  the  heat  pipe. 

9.4  START-UP  CONDITIONS 

Before  a  heat  pipe  can  reach  a  steady-state  condition,  the  energy  it  is  receives  as  heat 
input  must  heat  all  of  the  mass  that  comprises  the  heat  pipe  (fluid,  screen,  vessel  wall,  etc.) 
to  the  steady-state  temperatures.  The  USNA  heat  pipe  with  a  heat  input  of  400  W  would 
respond  as  shown  in  Figure  40.  Figure  40  compares  closely  to  Figure  37,  the  400W  start-up 
case  for  a  100%  charge.  The  experimental  case  reaches  a  slightly  higher  temperature  in  the 
activated  regions.  This  is  because  in  order  to  compensate  for  inactive  regions  which  are  at 
a  lower  temperature  the  active  regions  must  get  hotter  to  reject  the  input  heat.  However,  both 
graphs  reach  an  approximate  steady-state  condition  at  around  one  hour.  Therefore,  Figure 
40  represents  a  good  model. 


Transient  Start  Up  -  400W  ,  Horizontal 
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Figure  40.  Theoretical  Heat  Pipe  Start-Up 
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Another  aspect  of  startup  is  shown  in  Figure  41,  Figure  41  includes  a  plot  of  the 
heaters,  which  is  the  top  trace  in  this  graph.  It  is  interesting  to  note  that  when  heat  is  applied, 
the  heaters’  response  time  is  rapid.  The  heat  pipe  slowly  activates  more  regions,  and  when 
many  are  activated,  the  temperature  of  the  heater  drops  because  heat  is  being  more  effectively 
transferred  over  a  large  area. 
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Figure  41 .  Temperature  vs  Time  for  100%  Charge,  400-800W,  Horizontal 


9.5  ASYMMETRIC  CASES 

Two  cases  of  asymmetric  heating  will  be  discussed:  the  response  of  the  heat  pipe 
when  inclined,  and  the  response  when  horizontal  When  horizontal  and  at  a  low  heat  input, 
only  the  condenser  area  opposite  the  heater  section  activated,  as  seen  in  Figure  42.  However, 
as  the  heat  input  rises,  the  heat  pipe  condenser  region  begins  to  become  more  active  as  seen 
by  the  increasing  size  of  the  higher  temperature  section  of  the  condenser  surface.  Figures  43 
and  44  shows  the  heat  pipe’s  thermal  profile  when  the  heat  input  is  raised.  This  performance 
is  indicative  of  having  a  non-condensable  gas,  such  as  air,  in  the  heat  pipe.  If  the  USNA 
heat  pipe  performed  as  a  theoretically  perfect  heat  pipe  without  non-condensable  gases 
inside,  then  all  the  condenser  surface  in  the  asymmetric  (or  any)  case  would  reach  an 
isothermal  state. 
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Figure  42.  Asymmetric,  Horizontal,  200W 


Figure  43.  Asymmetric,  Horizontal,  400W 
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Figure  44,  Asymmetric,  Horizontal,  600W 


The  inclined  cases  did  not  perform  as  well  as  was  hoped.  The  condenser  side 
adjacent  to  the  heaters  never  flilly  activated,  and  the  more  developed  profile  remained  the 
“side-climbing”  phenomenon  that  was  observed  in  early  inclined  cases.  Figure  45  shows  the 
side  near  the  heater  section,  and  Figure  46  shows  the  larger  area’s  thermal  profile. 
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Figure  45  Asymmetric,  Inclined,  400W 
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Figure  46.  Asymmetric,  Inclined,  400W 


9.6  VARIABLE  CONDUCTANCE  HEAT  PIPES 

A  variable  conductance  heat  pipe  is  a  heat  pipe  that  operates  with  a  certain  amount 
of  non-condensable  gas  loading.  It  is  called  variable  conductance  due  to  the  effect  that  gas 
loading  has.  At  low  heat  inputs,  the  gas  blocks  off  a  large  portion  of  the  condenser  region, 
raising  the  temperature  for  the  part  that  can  be  utilized.  At  higher  powers,  less  of  the 
condenser  is  blocked  off  as  the  gas  is  compressed  more  against  the  end  of  the  heat  pipe. 

The  maximum  temperature  drop  across  the  wick  (including  liquid)  and  vessel  wall 
was  calculated  to  be  less  than  1 . 1  °C  for  a  heat  input  rate  of  500W.  Therefore,  we  cannot  use 
excess  lluid  to  explain  the  large  temperature  gradients.  Only  entrapped  air  in  the  heat  pipe 
could  be  causing  such  a  gradient  Observe  Figures  36-38,  and  note  that  the  end  temperatures 
all  remain  much,  much  lower  than  the  activated  regions.  These  are  areas  inside  the  heat  pipe 
that  are  covered  with  gas  and  the  water  vapor  could  not  push  the  gas  out  of  the  way  in  order 
to  have  an  active  condenser  section.  Therefore,  very  little  heat  is  transferred  in  that  region 
and  it  remains  much  colder 
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10.0  CONCLUSIONS 

Overall,  this  project  was  successful  in  determining  the  usefulness  of  flat  heat  pipes 
for  thermophotovoltaic  energy  conversion. 

10.1  GOALS  MET 

The  following  criteria  were  set  for  the  project: 

1.  Design  and  construct  the  heat  pipe  entirely  out  of  Monel  400. 

2.  Exhibit  two-dimensional  heat  transfer  behavior. 

3.  Feature  ultimate  dimensions  of  4  feet  long  (1 .22  m)  by  1  foot  wide  (30.48  cm)  by 
Vi  inch  thick  (1.27  cm). 

4.  Maintain  an  operating  temperature  of  up  to  212°F  (100°C). 

5 .  Possess  a  means  of  internal  support  for  structural  rigidity  without  sacrificing  two- 
dimensional  flow  capabilities. 

6.  Operate  without  the  presence  of  any  noncondensable  gases,  such  as  air,  inside  of 
the  heat  pipe  (gas  loading). 

The  project  had  the  following  experimental  and  theoretical  goals: 

1 .  Monitor  the  condenser  region  temperature  profile  as  closely  as  possible,  via 
thermocouples  and  infrared  video  camera. 

2.  Test  asymmetric  heating  conditions  and  differing  heat  inputs. 

3.  Test  performance  at  different  heat  inputs,  fluid  charges,  and  varying  inclination 
angles. 

4.  Modify  important  heat  pipe  equations  for  flat  geometry. 

5.  Model  experimental  and  theoretical  heat  pipe  performance  against  a  conventional 
fm  of  the  same  dimensions. 

6.  Analyze  experimental  data  to  determine  overall  performance  and  feasibility  of 
flat  heat  pipes  for  a  TPV  energy  conversion  system. 

The  goals  were  met  as  follows: 

•  A  flat  heat  pipe  of  the  set  dimensions  was  designed  entirely  of  Monel  400. 

•  The  flat  heat  pipe  was  successfully  constructed  using  pin  spacers  as  internal  supports 
and  operated  without  biasing  flow  direction. 

•  Using  water  as  a  working  fluid,  an  operating  temperature  in  excess  of  100°C  was 
achieved. 

•  The  heat  pipe’s  temperature  profile  was  monitored  and  recorded  closely. 
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•  Asymmetric  tests  were  conducted  and  clearly  demonstrated  that  a  heat  pipe  can 
maintain  an  isothermal  condenser  area  despite  a  nonuniform  heat  application. 

•  A  1 00%  charge  performed  the  best,  achieving  the  most  uniform  temperature  profile. 

•  Increasing  angle  of  orientation  from  horizontal  requires  a  higher  surface  temperature 
in  the  condenser  region  to  reject  the  same  amount  of  heat  the  vertical  orientation. 

•  The  horizontal  angle  of  orientation  performed  closest  to  a  theoretical  heat  pipe, 
achieving  the  most  near  isothermal  profile. 

•  Cylindrical  heat  pipe  equations  were  modified  to  flat  geometry  to  design  flat  heat 
pipes. 

•  Heat  pipes  were  modelled  against  fins  so  that  heat  pipe  advantages  could  be  better 
illustrated  as  an  isothermal  emitter  and  to  verify  experimental  operation  as  a  heat 
pipe 

•  The  heat  pipe  did  not  operate  completely  without  gas  loading.  However,  even  in  the 
presence  of  gas  loading,  it  is  possible  to  maintain  a  near  isothermal  condenser 
surface. 

•  It  has  been  demonstrated  that  a  flat  heat  pipe  can  maintain  a  large  isothermal  surface, 
with  either  a  uniform  (symmetric  case)  or  nonuniform  (asymmetric  case)  application 
of  heat. 

•  A  flat  heat  pipe  can  effectively  serve  as  an  emitter  for  a  thermophotovoltaic  energy 
conversion  system. 

10.2  RECOMMENDATIONS 

The  following  are  recommendations  based  on  the  data  gathered  by  this  project: 

•  Retest  heat  pipe  without  gas  loading. 

•  Test  higher  heat  fluxes  (near  capillary  limit)  to  verify  heat  pipe’s  full  heat  transport 
capabilities. 

•  Develop  transient  models  for  multidimensional  behaviors,  such  as  large  temperature 
gradient  in  inclined  cases  and  temperature  oscillations. 
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APPENDIX  1  -  CAPILLARY  LIMIT  DERIVATIONS 

Flat  Heat  Pipe 


70 


Capillary  Limit  (neglect  boiling  incipience)  for  horizontal  case 


AP  =  P  -  P  = 

^■'^max  ■‘v  ■'Z, 


2g 


applying  the  Clapeyron  relationship; 


dP\ 


dT) 


sat 


since  Vf  is  small,  Vfg  =  Vg 


therefore 


dP  h 

dT  rVj  RT^ 


integrating  from  P^,  Tf  to  Pl,  Tl 


In 


fAl 

1 

R 

T 

n 


recall 


P  =  P. 


sat 


2(j 
+  — 
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Pressure  in  Flat  Heat  Pipe 


"^sat 

Vapor 

Wick 

~P,at 

Liquid  -  Vapor  Path 


AP  =AP  +APr  . 

max  ^  vapor  liquid 


but  the  path  is  some  averaged  path  estimate  of  averaged  path: 


-  L  L 

^  =  -:r+4  +  vf 

2  “  2 


for  no  adiabatic  section,  as  in  USNA  fhp 


dP  ^  d'u 

_  —  =  G= 
dx  dy 


9yL 

2// 


+  Ay  B 


Roimdarv  Conditions 
at  y=0,  u=0 


K  ou 

at  J  =  0 

2  ay 


Applying  Boundary  Conditions, 


then, 


K 


fmdP;^  = 


12// 
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so, 


dx 


h; 


L 


let  RCy-  = 


film  Reynold’s  number 


12 

Re^ 


Pressure  Loss  Through  Vapor  assumino  laminar  flow 


^P.. 


24  L 

Re^ 


fn  =  P.Kd^v  = 


where  -  Wh^  (Vapor  Flow  Area) 


^  ^  \2q,Lp^ 

^  ■  pXA  ■  PXAh,, 

where 

hv  =  6v 
So, 


AR, 


PXAhfs 


A  A 
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Vapor  Flow  (Find  effective 


y 


z 


Lt 


.56 


V 


Lc 


Assume  Re^  = 


p  5V 

rv  V  gy 

Av 


let 


dz) 


\lq,p,  llmp^ 


PA 


pA  a 


for  mass  flow  in  z  direction, 


^  <7. 

m.  =  —^,0<z<  L 
-Jh  ^ 

^fg 


' 


= 


Li'j'  z 


4 -(4 +  4) 


f  A 


yK/ 


,  Lg  +  ^  z  Lj- 
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dP  z  (Lj-  -  z) 

—  =-B  — +  1+ — — — 

4  4 -(4 +  4) 


where 


integrating  both  sides, 


p,(Lj)-pm--b<  jj-  +Wi;*‘-+  4^ 

L  ^-10  L  ^  J/  + 


this  simplifies  to: 


PALr)-PM  = 


( L 


/^v4  44  ^  2  2 


■+  + 


which  is  the  pressure  loss  in  the  vapor  space  where: 


4  4l  - 

—  +  4  +  —  =  L 

2  "27 


Flow  through  Porous  Media 


For  Flat  Plate  (non-porous)  channel: 


AP, 


llmLfi 

p5^A 


A  =  W6 


recall, 


therefore 


12  12// 
Re^  "  pV5 

AP  \lpV 
L  ~  5^ 
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and 

dP  \2fiV 
dx  5~ 

For  porous  media, 


dx  PlA^^K 

K  =  permeability  (m^) 
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simplifies  to  approximately 


set  this  equal  to  the  earlier  result, 

dP_ 

dx  PiAJi 


solving  for  K, 


p,A,K 


PKifiT) 


however,  using  permeability  equation  for  cylindrical  porous  media  and  approximating 
with  Dh,  hydraulic  radius  is  a  valid  approximation: 


to  find  liquid  pressure  drop  in  heat  pipes,  use  this  K  along  with  the  following  equation: 

dP  _ 

dx  p^AJPi 
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Capillary  Limit  -One  Side  heated 


Steady  state  limit 


max  capillary  pumping  head  (neglect  gravity) 
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AP 


cap 


dP 

+  —~dx 
dx 


AP 


cap 


12//>  - 


2(7 


/2 


L\ 


'fg 


12Av 


■  + 


Ml 


[p^A^S:  p,A,K\ 


therefore,  capillary  limit,  one  side  heated  is: 


)  _  2“"^*  I 

12a  ^  Ml 

'c°Phm„  ^effL 

[p^A^Sp  p, A, k\ 
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and 


Capillary  Limit.  Both  Sides  Heated 

AP,  =  f,m 

APi.1  -  f L\^L\ ’ 

APz,2  ~  f Ll^Ll 
APj.  =  fjm  =  f^m\ 

APr  =  fr^  =  fiei^Lx  +  ^Li) 


letting  APl,  =  APl2 


where 


and 


therefore 


A^r  -  AP,  =  /„AP, 


M  1  ^ 


+ 


Ll^ 


f-^L 


fix  ~  f LI  ~ 


Pi 


Pl^l^ 


L 


AP. 


f/v  +  P 

=  ihL 

12Pv_^  Pi 

O  ' 


p,A,S;  2p,A,K 


and,  Capillary  Limit  for  2  sides  heated: 


12//,  ^  Ml 

^effL  '^Pl^lK 
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Appendix  2  -  BOILING  LIMIT  DERIVATION 


Liquid 


1 


dP  = 


2g 


Fb  =  effective  bubble  radius 


Clapeyron  Equation; 


'^fg  (small)  =  Vg 

v^-Vb=  (RTb)/Pb,  Pb  =  P 

_  V 

dT)  "  RT^ 

sat 


I 


dT 

j2 


^satj 


recall, 


so 


where 


2r7 


2(j 

Psatf 


h 


'fg 

R 


Psat,l  Pb 
PsaljPb 


T  X  T  X  T 

^sat,l  ~^b  ^sat 


Therefore,  wall  temperature  difference  required  for  nucleation 


RlcjT 


sat  A 


nt  ' 


Z  -  Z 
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Screen 


Vapor  Area 

T  aat- 


Liquid 


T  *iah  1 


Bubble 
Pb.  Tb 


T„  =  Tb 


Note:  rh«ra, 


Vessel 

material 

(wall) 


Applied  to  heat  pipe  wick 

Boiling  Onset  in  Heat  Pipe 
AP  from  vessel  wall  to  vapor  area: 

Pt  -  c  =  (p,  -  -  />„,,) 


P,-P^.-2a[j;--i] 


Using  Clapeyron  Criteria  for  boiling 


IcrRr, 

AZ.:,  =  L.  -  = 


^sat^fg 


Heat  into  heat  pipe  is  transferred  by  conduction, 

kAAT 

Qcond  - 


therefore,  boiling  limit  is  reached  when 


a  = 


where  U,  the  overall  heat  transfer  coefficient  is  found  by: 


U  = 


where  C  is  vessel  thickness  and  U  is  wick  thickness 


therefore. 


2 

sat 


P  h 


1 


vt' 


'  + 


APPENDIX  3  -  SONIC  LIMIT  DERIVATION 


at  steady-state: 

Q  =  rhhj-^ 


where  th  =  pA  V 

also  max  velocity  (sonic)  for  constant  area  duct  is: 


local 


Qsonic  = 


substituting. 
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APPENDIX  4  -  DEFLECTION  DERIVATIONS  AND  SPREADSHEET 


Pin  to  Pin  Case 

iP  =  14.7s 


Me! 


c 

'VA 


W=  14,7  lb/  in  for  unit 
inch  width 


Expected  deformation 


. A 

R. 


s  =  spacer  distance  from  wall 


y  A 


~  (symmetry) 


A 


X* 


d~v 

Equation  1;  EI-^^=  R^x  -  7 35x 


Equation  2;  El x^  -  ^ x^  +  C. 
ax 


Equation  3  :  EIv  =  +  C^x  +  Cj 


Boundary  Conditions 
at  X  =  0,  V  =  0.  Therefore  Cj  =  0 
at  X  =  s/2,  dv/dx  =0 


C,  =  -.61255' 

knowing  that  is  at  x  =  s/2  and  solving  equation  3; 


V 


max 


-.1914055^ 

El 
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Sidewall  to  Pin  Case 

iP  =  14 . 7s 


W=  14.7  lb/  in  for  unit 
inch  width 


VVVVl  ^ 

.  Expected  deformation 

m----'- . . - . - . > 

A  s  =  spacer  distance  from  wall 


Equation  1  M  4-  V^x  - 

ax 


dv  K  ,  Wx^ 

Equation  2  El  —  =  M^x  +  —  x  - +  C 

<fc  "  2  6 


K  2  K  ^  ^4 

Equation  3  EIv  = - x  +  — x'’  -  — x  +  C,x  +  C, 

2  6  24  '  ^ 


Boundary  Conditions 
at  X  =  0,  V  =  0,  dv/dx  =  0 
therefore  C,  =  0  =  C2 

To  find  Me  and  use  equations  1  and  3  and  at  A,  M  =  0  and  v  =  0 


Mc  =  - 1.8375^ 


Vc  =  9.18755 
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solving  for  max  deflection  (dv/dx  =  0)  using  equation  2 

at  X  =  .57846483s 


substituting  into  equation  3 


V 


max 


.22785/-.30747.!>^ 

El 
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Rigid  Pin  to  Rigid  Pin  Case 
(welded  both  sides) 


Using  Ferdindand  Singer’s  Strength  of-Materials  (1962); 


Ws^ 

12 


And 


_ 

“  384£/ 
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Materials 

psipsi 

EG 

psi  psi 

Yield  T  Yield  C 

psi  psi 

UItT  UltC 

Poisson's 

Monel 

R-400 

2.6E+07  9.4E+06 

3.5E+04 

8.0E+04 

0.295 

R-405 

2.6E+07  9.4E+06 

7.4E+04 

8.3E+04 

0.32 

K-500 

2.6E+07  9.4E+06 

1.2E+05 

1.6E+05 

0.32 

SS  304 

2.80E+07 

5.50E+04 

1.00E+05 

0.3 

1040 

3.00E+07 

5.43E+04 

8.55E+04 

0.3 

for  unit 
inch  width 

Desired  0.0625 

deflection  = 

Spacing  1040 

Material  Thickness 

0.035  0.04  0.045  0.05  0.055  0.06  0.065 

0.1  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.2  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.3  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.4  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.5  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.6  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.7  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.8  -0.0001  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

0.9  -0.0001  -0.0001  -0.0000  -0.0000  -0.0000  -0.0000  -0.0000 

1  -0.0001  -0.0001  -0.0001  -0.0000  -0.0000  -0.0000  -0.0000 

1.1  -0.0002  -0.0001  -0.0001  -0.0001  -0.0001  -0.0000  -0.0000 

1.2  -0.0003  -0.0002  -0.0001  -0.0001  '  -0.0001  -0.0001  -0.0000 

1.3  -0.0004  -0.0003  -0.0002  -0.0001  -0.0001  -0.0001  -0.0001 

1.4  -0.0005  -0.0003  -0.0002  -0.0002  -0.0001  -0.0001  -0.0001 

1.5  -0.0007  -0.0005  -0.0003  -0.0002  -0.0002  -0.0001  -0.0001 

1.6  -0.0009  -0.0006  -0.0004  -0.0003  -0.0002  -0.0002  -0.0001 

1.7  -0.0011  -0.0007  -0.0005  -0.0004  -0.0003  -0.0002  -0.0002 

1.8  -0.0014  -0.0009  -0.0007  -0.0005  -0.0004  -0.0003  -0.0002 

1.9  -0.0017  -0.0012  -0.0008  -0.0006  -0.0004  -0.0003  -0.0003 

2  -0.0021  -0.0014  -0.0010  -0.0007  -0.0006  -0.0004  -0.0003 

2.1  -0.0026  -0.0017  -0.0012  -0.0009  -0.0007  -0.0005  -0.0004 

2.2  -0.0031  -0.0021  -0.0015  -0.0011  -0.0008  -0.0006  -0.0005 

2.3  -0.0037  -0.0025  -0.0018  -0.0013  -0.0010  -0.0007  -0.0006 

2.4  -0.0044  -0.0030  -0.0021  -0.0015  -0.0011  -0.0009  -0.0007 

2.5  -0.0052  -0.0035  -0.0025  -0.0018  -0.0013  -0.0010  -0.0008 

2.6  -0.0061  -0.0041  -0.0029  -0.0021  -0.0016  -0.0012  -0.0010 

2.7  -0.0071  -0.0048  -0.0033  -0.0024  -0.0018  -0.0014  -0.0011 

2.8  -0.0082  -0.0055  -0.0039  -0.0028  -0.0021  -0.0016  -0.0013 

2.9  -0.0095  -0.0063  -0.0045  -0.0032  -0.0024  -0.0019  -0.0015 

3  -0.0108  -0.0073  -0.0051  -0.0037  -0.0028  -0.0022  -0.0017 

3.1  -0.0124  -0.0083  -0.0058  -0.0042  -0.0032  -0.0025  -0.0019 
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APPENDIX  5  -  MATHCAD  MODELLING  OF  FIN  VS  HEAT  PIPE 


Fin  Analysis  -  convection  only  -  steady  state 

k:=150 


Ax:=—l 

12 


:=1.5 

P:=2  +  -L 

12 

._/h-Pr 
♦  ~  - 

T 

room 

\k-Axi 

:=70 


Q.  ;=  100-3. 41 2j 


Qj 

Qwatt.  :=-  ^ 


J  3.412 


TO',  x)  := 


T  -H 
^  room^ 


j 


k-Ax-ni/ 


•exp(-m-x) 


-32 


Distance  from  Evap,  ft 


95 


Tsat,  isothermal  steady  state  heat  pipe  analysis  with  radiation  and  convection. 

0  :=0.171810*^  E  :=0.95  As  :=6 

Qconv(T33)  :=h  As.(T33-T,„„^k)  Q'-^d(T33)  :=(a  e -As ).(t 33^- T 

T33  :=460 

T3i(j):=root(Q^-Qconv(T33)-Qrad(T33).T33)  Ts.  :=[(T 3i(j)- 460)  -  32].l 

_Qconv(T3j(j))  ^  _Qrad(T3i(j)) 

I4T2  1412 


0  100  200  300  400  500  600  700  800  900  1000 

Qwatt. 

} 

Power  In,  W 


convection 

radiation 
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Steady  state  heat  pipe,  combined  convection  and  radiation  modes,  variation  of  heff  with  angle  of 
inclination.  Q  in  =  400W 

n  :=0..  18  r4.124 

u  .  X  -  4.120 

theta  :=(n)o 


heff  :=  3.781 


Qconvl(T,,.n)  -(heff 

3.2S- 


.c.  Ts2  anc 


Qradl(T33):=(a.E-As)T33^-T,,,^k' 


Tss  -460 


T  s2(  n  )  •=  root  ^Q4  —  Qconv  1  [T  ,  n j  -  Qrad  1 

T,2_a„g„:=[(T,2(n)-460j-32]i 
Q,  =  1.36510^ 

4 


Variation  of  Thp  with  angle  for  400 W  Qin 


10  20  30  40  50  60  70  80  90 


Angle,  degrees 


;_ang  = 


57.023 

57.045 

57.111 

57.222 

57.378 

57.587 

57.843 

58.162 

58.544 

59.008 

59.556 


61.005 

61.984 

63.214 


Transient  lumped  parameter  heat  pipe  fin  model-  Monel  and  water  masses  are 
combined  with  appropriate  specific  heat.  The  heat  Pipe  temperature  is  assumed  to 
be  isothermal. 


0=m’C  — Trad(t)  -  Q -t- h'A-(Trad(t)  -  TreO  -i- A-a-s-(^Trad(t)'^  -  Tref^) 
dt 


m  :=  3.85  c  :=  1  h  ;=  1.23 


1.365  10^  As  =6 


Thpl  :=  530,535..  620 


rThpl 


t2{Thpl)  ;=  m  e- 


J530 


1 

-  h  As'(Trad  -  530)  -  As-a-E-iTrad"*  -  530'*) 


dTrad 


ThpC(Thpl)  :=  (Thpl  -  460  -  32)--- 


'fhp  I  ) 
rime,  hr 
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Fin  Analysis  -  convection  only  (assume  radiation  losses  are  close  to  convection  losses)  see  other 
analysis  -  steady  state,  so  h  Is  approximately  3, 


j:=l,.  10 


h:=3 


P:=2+- 


k  :=  150 
1 

12 


Ax  :=  — •  1 
12 


m  := 


h'P 


0.5 


\k-Ax 
Thpipe(x)  :=  68. 


T  :=  70 
^  room 


100  3.412  j  Q, 


J  3.412 


T(j.x)  := 


-32 


cm(x)  := - •  100 

3-28 


cm(x) 

Distance  from  Evap,  cm 


Regular  Fin,  400 W 
Heat  Pipe  Fm,400W 
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